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Q Abstract: The study of finite J effects for the hght-cone AdSs x superstring 

by means of the Thermodynamic Bethe Ansatz requires an understanding of a com- 
,—1 panion 2d theory which we call the mirror model. It is obtained from the original 

string model by the double Wick rotation. The S-matrices describing the scattering 
p I of physical excitations in the string and mirror models are related to each other by 

an analytic continuation. We show that the unitarity requirement for the mirror S- 
matrix fixes the S-matrices of both theories essentially uniquely. The resulting string 
CO S-matrix S{zi,Z2) satisfies the generalized unitarity condition and, up to a scalar 
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J> solutions, however, exhibit various signs of pathological behavior. In particular, they 

^ might receive a finite J correction to their energy which is complex, or the energy 

correction might exceed corrections arising due to finite J modifications of the Bethe 
equations thus making the asymptotic Bethe ansatz inapplicable. 
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1. Introduction and summary 

The conjectured duahty between the maximally supersymmetric Yang-Mills theory 
in four dimensions and type IIB superstring in the AdSs x background [1] is the 
subject of active research. Integrability emerging on both sides of the gauge/string 
correspondence [2, 3] proved to be an indispensable tool in matching the spectra of 
gauge and string theories. Namely, it was shown that the problem of determining 
the spectra in the large volume (charge) limit, can be reduced to the problem of 
solving a set of algebraic (Bethe) equations. The corresponding Bethe equations are 
based on the knowledge of the S-matrix which describes the scattering of world-sheet 
excitations of the gauge-fixed string sigma-model, or alternatively, the excitations of 
a certain spin chain in the dual gauge theory [4]- [9]. 

Remarkably, the S-matrix is severely restricted by the requirement of invariance 
under the global symmetry of the model, the centrally extended psu(2|2) © psu(2|2) 
superalgebra [9, 10]. The invariance condition fixes its matrix form almost uniquely 
up to an overall phase [9, 11]. The constraints on the overall phase were derived 
in [12] by demanding the S-matrix to satisfy crossing symmetry. Recently, a physi- 
cally relevant solution to the crossing relation, which interpolates between the weak 
(gauge) and strong (string) coupling regimes was conjectured [13, 14], building on 
the previous work [15]- [18]. This solution successfully passed a number of non-trivial 
tests [19]- [34]. 

So far the main focus of research was on determining the spectrum of string 
theory in the limiting case in which at least one of the global charges carried by a 
string state (and by the corresponding gauge theory operator) is large. Our ultimate 
goal, however, is to understand how the energies of string states (the conformal 
dimensions of dual gauge theory operators) depend on the coupling constant for 
finite values of all the other global symmetry charges. Although the conjectured 
S-matrix [13, 14] provides an important starting point in addressing this issue, by 
now there is firm evidence that the corresponding Bethe equations [8] fail to correctly 
reproduce the finite-size effects, neither in string [35] nor in gauge theory [36]. Indeed, 
already in the semi-classical string theory deviations from the controllable exact 
spectrum arise which are exponentially small in the effective string length playing 
the role of a large symmetry charge. One of the reasons behind this is that the 
interactions on the world-sheet are not ultra-local, typically the scattered states are 
solitons of finite size [37, 38]. Also, as is common to many field-theoretic models, 
vacuum polarization effects smear bare point-like interactions and lead to exponential 
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corrections to the energy levels in the large- volume limit [39]. Complementary, in 
the spin chain description of the dual gauge theory, the obstruction to the validity of 
the Bethe ansatz comes from the wrapping effects. Hence, the Bethe ansatz for the 
planar AdS/CFT system in its present form is merely of asymptotic type. Obviously, 
solving a quantum sigma-model in finite volume is much harder. As an illustrative 
example, we mention the Sinh-Gordon model for which the equations describing the 
finite-size spectrum have been recently obtained in [40]. 

Determination of finite-size effects in the context of integrable models is a wide 
area of research. Basically, there are three different but related ways in which this 
problem could be addressed: the thermodynamic Bethe ansatz approach (TBA) 
[41], nonlinear integral equations (NLIE) [42] and functional relations for commuting 
transfer-matrices [43] . 

The aim of the present paper is to investigate the structure of the string S-matrix 
which is the necessary step for constructing the TBA equations for the quantum 
string sigma-model on AdSs x background. This would eventually allow one to 
describe the finite-size spectrum of the corresponding model. 

The TBA approach was initially developed for studying thermodynamic proper- 
ties of non-relativistic quantum mechanics in one dimension [44] and further applied 
to the computation of the ground state energy in integrable relativistic field theories 
in finite volume [41]. The method also was later extended to account for energies of 
excited states [45] (see also [46]). 

Implementation of the TBA approach consists of several steps. The primary 
goal is to obtain an expression for the ground state energy of a Lorentzian theory 
compactified on a circle of circumference L and at zero temperature. The starting 
point is the Euclidean extension of the original theory, put on a torus generated 
by two orthogonal circles of circumferences L and R. The partition function of 
this theory can be viewed as originating from two different Lorentzian theories: the 
original one, which lives on a circle of length L at temperature T and has the 
Hamiltonian if, or the mirror theory which is defined on a circle of length R = 1/T 
at temperature T = 1/L and has the Hamiltonian H. For Lorentz-invariant theories, 
the original and the mirror Hamiltonians are the same. However, in general and in 
particular for the case of interest here, the two theories need not be the same. Taking 
the thermodynamic limit R oo one ends up with the mirror theory on a line and 
at finite temperature, for which the exact (mirror) Bethe equations can be written. 
Thus, computation of finite-size effects in the original theory translates into the 
problem of solving the infinite volume mirror theory at finite temperature. 

Although taking the thermodynamic limit simplifies the system, a serious com- 
plication arises due to the fact that the mirror theory could have bound states which 
manifest themselves as poles of the two-particle mirror S-matrix. Thus, the com- 
plete spectrum would consist of particles and their bound states; the latter should 
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be thought of as new asymptotic particles. Having identified the spectrum, one has 
to determine the S-matrix which scatters all asymptotic particles. It is this S-matrix 
which should be used to formulate the system of TEA equations. 

As is clear from the discussion above, the mirror theory plays a crucial role in 
the TEA approach. In this paper we will analyze the mirror theory in some detail. 
First, we will explain its relation to the original theory. Indeed, given that the model 
in question is not Lorentz invariant, the mirror and the original Hamiltonians are not 
the same. However, since the dispersion relation and the S-matrix can be deduced 
from the 2-point and 4-point correlation functions on the world-sheet, and since the 
correlation functions in the mirror theory are inherited from the original model by 
performing a double Wick rotation, it follows that the mirror dispersion relation 
and the mirror S-matrix can be obtained from the original ones by the double Wick 
rotation. Here we will meet an important subtlety. To explain it, we first have to 
recall the basic properties of the string S-matrix. 

As was shown in [47], the psu(2|2) ©psu(2|2)-invariant S-matrix S{pi,p2), which 
depends on real momenta pi and p2 of scattering particles, obeys 

• the Yang-Eaxter equation 

»523»5i3iSi2 = 5121S131S23 

• the unitarity condition 

Sl2{Pl,P2)S2l{p2,Pl) = I 

• the physical unitarity condition 

Sl2iPl,P2)Sl2ipi,P2) = I 

• the requirement of crossing symmetry 

^rl5^2(pl,P2)^l5i2(-pi,P2) =1, 

where ^ is the charge conjugation matrix. 

The first three properties naturally follow from the consistency conditions of the 
associated Zamolodchikov-Faddeev (ZF) algebra [48, 49], while the last one refiects 
the fact that the particle-to-anti-particle transformation is an automorphism of the 
ZF algebra [47] . The unitarity and physical unitarity conditions imply the following 
property 

S2l{p2,Pl) = Sl2ipi,P2) ■ 

One should bear in mind that the S-matrix is defined up to unitary equivalence 
only: unitary transformations (depending on the particle momentum) of a basis of 
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one-particle states correspond to unitary transformations of the scattering matrix 
without spoiling any of the properties listed above. 

The mirror S-matrix S{pi,p2) is obtained from S{pi,p2) by the double Wick 
rotation. The above-mentioned subtlety lies in the fact that only for a very special 
choice of the one-particle basis the corresponding mirror S-matrix remains unitary. 
As we will show, this problem can be naturally attributed to the properties of the 
double Wick rotation for fermionic variables. Upon the basis is properly chosen to 
guarantee unitarity of the mirror S-matrix, the only freedom in the matrix structure of 
S{pi,p2) reduces to constant, i.e. momentum-independent, unitary transformations^. 

There is another interesting explanation of the interrelation between the original 
theory and its mirror. As was shown in [12], the dispersion relation between the 
energy and momentum of a single particle can be naturally uniformized in terms of 
a complex variable z living on a torus with real and imaginary periods equal to 2lji 
and 2uj2, respectively. Since z plays the role of the generalized rapidity variable, it 
is quite natural to think about the S-matrix as the function S{zi, Z2), which for real 
values of the generalized rapidity variables coincides with S{pi,p2). In other words, 
the S-matrix admits an analytic continuation to the complex values of momenta. It 
appears that the unitary momentum-dependent freedom in the matrix structure of 
the S-matrix gets fixed if we require the analytic continuation to be compatible with 
the requirement of 

• generalized unitarity 

Si2{zl,Z;)[S,2{Zl,Z2)y =1, 

which can be thought of as the physical unitarity condition extended to the general- 
ized rapidity torus. The unitarity and the generalized unitarity further imply 

S21{Z2,ZI) = [Si2{Zi,Z2)y 

In fact, the last equation is equivalent to the standard requirement of hermitian 
analyticity for an S-matrix in two-dimensional relativistic quantum field theories. 

Thus, the S-matrix which admits the analytic continuation to the generalized 
rapidity torus compatible with the requirement of hermitian analyticity is essentially 
unique. Of course, it satisfies all the other properties listed above, including crossing 
symmetry. As we will show, the mirror S-matrix is obtained from S{zi, Z2) considered 
for real values of Zi, Z2 by shifting zi, Z2 by quarter of the imaginary period 

S{ZI,Z2) =S{zi + f,Z2+f) . 

^Of course, there is always a freedom of multiplying the S-matrix by an overall (momentum- 
dependent) phase. 
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There is a close analogy with what happens in relativistic models. In the latter 
case the physical region is defined as the strip < ImO < vr, where 6 = 62 — 6i 
is the rapidity variable. A passage to the mirror theory corresponds to the shift 
Ok ^ Ok + if, i.e. to the shift by the quarter of imaginary period^. Of course, for 
relativistic models, due to Lorentz invariance, the S-matrix depends on the difference 
of rapidities and, therefore, it remains unchanged under the double Wick rotation 
transformation. Also, in our present case the notion of the physical region is not 
obvious and its identification requires further analysis of the analytic properties of 
the string S-matrix. 

Having identified the mirror S-matrix, we can investigate the question about the 
bound states. We first discuss the Bethe equations for the gauge-fixed string theory 
where the existence of the BPS bound states is known [50]. No non-BPS bound states 
exist, according to [50]-[52]. We find out, however, that the number of solutions of the 
BPS bound state equations depends on the choice of the physical region of the model, 
and for a given value of the bound state momentum there could be 1, 2 or 2^^"^ M- 
particle bound states sharing the same set of global conserved charges. It is unclear 
to us whether this indicates that the actual physical region is the one that contains 
only a single M-particle bound state or it hints on a hidden symmetry of the model 
responsible for the degeneracy of the spectrum. These solutions behave, however, 
differently for very large but finite values of L; most of them exhibit various signs of 
pathological behavior. In particular, they might have complex finite L correction to 
the energy, or the correction would exceed the correction due to finite L modifications 
of the Bethe equations thus making the asymptotic Bethe ansatz inapplicable. In 
the weak coupling limit, i.e. in perturbative gauge theory, and for small enough 
values of the bound state momentum only one solution reduces to the well-known 
Bethe string solution of the Heisenberg spin chain. It is also the only solution that 
behaves reasonably well for finite values of L. Therefore, it is tempting to identify 
the physical region of the string model as the one that contains this solution only. 

By analyzing the Bethe equations for the mirror theory, we show that bound 
states exist and that they can be regarded as "mirror refiections" of the BPS bound 
states in the original theory. No other bound states exist, in agreement with the 
results by [52]. Given the knowledge of bound states, the next step would be to 
construct the S-matrix which describes scattering of all asymptotic particles including 
the ones which correspond to bound states. In principle, such an S-matrix can be 
obtained by the fusion procedure [53, 54] applied to the "fundamental" S-matrix we 
advocate here. This is the bootstrap program whose discussion we will postpone for 
the future. 

The paper is organized as follows. The next section contains the discussion of 
the double Wick rotation, the mirror dispersion relation and the mirror magnon. In 

^The shift of 6 by the half-period corresponds to the crossing transformation. 
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section 3 we discuss the supersymmetry algebra and the construction of the mirror S- 
matrix. In section 4 we analyze the double Wick rotation on the generalized rapidity 
torus as well as various possible definitions for the physical region. In section 5 
the properties of the string S-matrix defined on the generalized rapidity torus are 
discussed. We also prove here the unitarity of the scalar factor in the mirror theory. 
In section 6 we present various versions of the Bethe equations in the original and 
mirror theory pointing out that the Bethe equations based on the 3u(2|2) ©su(2|2)- 
invariant string S-matrix should be modified in the odd winding number sector since 
for this case the fermions of the gauge-fixed string sigma model are anti-periodic. 
Sections 7 and 8 contain an analysis of the bound states of the AdSs x S^ gauge-fixed 
model and its mirror theory. Section 9 consists of several appendices. 



2. Generalities 

In this section we discuss how the vacuum energy of a two-dimensional field theory 
on a circle can be found by considering the Thermodynamic Bethe Ansatz for a 
mirror model obtained from the field theory by a double Wick rotation. We follow 
the approach developed in [41]. 

2.1 Double Wick rotation and mirror Hamiltonian 

Consider any two-dimensional field-theoretic model defined on a circle of circumfer- 
ence L. Let 

H= I danip,x,x') (2.1) 
Jo 

be the Hamiltonian of the model, where p and x are canonical momenta and coor- 
dinates. They may also include fermions but in this section we confine ourselves to 
bosonic fields only. We will refer to the action corresponding to the Hamiltonian H 
as to the Minkowski action, however it does not have to be relativistic invariant. 

We want to compute the partition function of the model defined as follows 

Z{R, L) ^ 5^(^„|e-^^|^„) = ^ e-^"^ , (2.2) 

n n 

where \il)n) is the complete set of eigenstates of H. By using the standard path 
integral representation [55], we get 

Z{R,L) = J VpVxe^o'^'i^fo Mip^-n) ^ (2.3) 

where the integration is taken over x and p periodic in both r and a. Formula 
(2.3) shows that — dr da{ipx — Ti.) can be understood as the Euclidean action 
written in the first-order formalism. Indeed, integrating over p in the usual first-order 
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action dr da{px — Ti), we get the Minkowski-type action, and the Euchdean 
action is obtained from it by replacing i; — > which is equivalent to the Wick 
rotation r —it. 

Let us now take the Euclidean action, and replace x' —ix' or, equivalently, do 
the Wick rotation of the cr-coordinate a — > ia. As a result we get the action where 
(7 can be considered as the new time coordinate. Let H be the Hamiltonian with 
respect to a 

H= I dTn{p,x,x), (2.4) 
Jo 

where p are canonical momenta of the coordinates x with respect to a. 

We will refer to the model with the Hamiltonian H as to the mirror theory. If 
the original model is not Lorentz-invariant then the mirror Hamiltonian is not equal 
to the original one, and the Hamiltonians H and H describe different Minkowski 
theories. 

The partition function of the mirror model is given by 

Z{R, L) ^ J](V^Je-^^|V^„) = ^"^"^ ' (2.5) 

n n 

where j^'n) is the complete set of eigenstates of H. Again, by using the path integral 
representation, we obtain 

Z{R, L) = J VpVx e^o"" dr da(ipx'-n) _ ^2.6) 

Finally, integrating over p, we get the same Euclidean action and, therefore, we 
conclude that 

Z{R,L) = Z{R,L) . (2.7) 

Now, if we take the hmit it! — > 00, then log Z{R, L) ~ —RE{L), where E{L) is the 
ground state energy. On the other hand, \ogZ{R, L) ^ —RLf{L), where f{L) is the 
bulk free energy of the system at temperature T — 1/L with a considered as the 
time variable. This leads to the relation 

E{L) = Lf{L) . (2.8) 

To find the free energy we can use the thermodynamic Bethe ansatz because i? >> 1. 
This requires, however, the knowledge of the S-matrix and the asymptotic Bethe 
equations for the mirror system with the Hamiltonian H. Although the light-cone 
gauge-fixed string theory on AdS^ x is not Lorentz invariant, H ^ H, it is still 
natural to expect that there is a close relation between the two systems because their 
Euclidean versions coincide. 
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A potential problem with the proof that Z{R, L) = Z{R, L) is that the integra- 
tion over p and p produces additional measure factors which may be nontrivial. The 
contribution of such a factor is however local, and one usually does not have to take 
it into account. We will assume throughout the paper that this would not cause any 
problem. 

2.2 Mirror dispersion relation 

The dispersion relation in any quantum field theory can be found by analyzing the 
pole structure of the corresponding two-point correlation function. Since the corre- 
lation function can be computed in Euclidean space, both dispersion relations in the 
original theory with H and in the mirror one with H are obtained from the following 
expression 

i^i + Vsin^f + 1, (2.9) 

which appears in the pole of the 2-point correlation function. Here and in what 
follows we consider the light-cone gauge-fixed string theory on AdSs x which has 
the Euclidean dispersion relation (2.9) in the decompactification limit L = P+ — > oo 
[5, 6, 9, 10]. The parameter g is the string tension, and is related to the 't Hooft 
coupling A of the dual gauge theory as g = 

Then the dispersion relation in the original theory follows from the analytic 
continuation (see also [39]) 

He^-iH, pe-^p ^ i72_l + Vsin2|, (2.10) 

and the mirror one from 

He^P, PE-^iH H = 2 arcsinh (^^v^l +p2j . (2.11) 

Comparing these formulae, we see that p and p are related by the following analytic 
continuation 



p ^ 2i arcsinh (^^V^TTF) , H = ^ 1 + Ag"^ sin^ ^ ^ ip . (2.12) 

We note that the plane-wave type limit corresponds to taking g ^ oo with p fixed, 
in which case we get the standard relativistic dispersion relation 

^pw = ^A/r+F- (2.13) 

The expression above suggests that in this limit it is natural to rescale H by 1/g 
or, equivalently, to rescale t = ia hj g. This also indicates that the semi-classical 
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limit in the mirror theory should correspond to (7 — > cxd with p/ g fixed, so that the 
dispersion relation acquires the form 

^sc = 2arcsinh (^^^ . (2.14) 

We will show in the next subsection that the mirror theory admits a one-soliton 
solution whose energy exactly reproduces eq.(2.14). 

In what follows we need to know how the parameters introduced in [5] are 
expressed through p. By using formulae (2.12), we find 




^ (P)^-A\l^ + V^. T l\{p-t) (2.15) 
and, as a consequence, 

ix~ — ix~^ —^-(l + ip) . 
9 

Note that these relations are well-defined for real p, but one should use them with cau- 
tion for complex values of p. In section 4 we introduce a more convenient parametriza- 
tion of the physical quantities in terms of a complex rapidity variable z living on a 
torus [12]. In this parametrization the analytic continuation would simply correspond 
to the shift of z by the quarter of the imaginary period of the torus. 

2.3 Mirror magnon 

In this section we will derive the dispersion relation for the "giant magnon" in the 
semi-classical mirror theory. This will provide further evidence for the validity of the 
proposed dispersion relation (2.14). 

Consider the classical string sigma-model on AdSs x and fix the generalized 
uniform light-cone gauge as in [56, 57]. The gauge choice depends continuously on a 
parameter a with the range < a < 1. The gauge-fixed Lagrangian in the generalized 
a-gauge can be obtained either from the corresponding Hamiltonian [59, 57] by using 
the canonical formalism or by T-dualizing the action in the direction canonically 
conjugate to the light-cone momentum P+ [26]. Its explicit form in terms of the 
world-sheet fields is given in appendix 9.1. To keep the discussion simple, in what 
follows we will restrict our analysis to the a = 1 gauge^. 

We are interested in finding a soliton solution in the mirror theory, which is 
obtained from the original theory via the double Wick rotation with further exchange 
of the time and spacial directions 

a = —ir , f = ia , (2.16) 

^Recall that unlike to the case of finite P+ the dispersion relation of the giant magnon in the 
infinite volume limit P+ = 00 was shown to be gauge independent [38]. 
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where a, r are the variables parametrizing the world-sheet of the original theory. 

Recall that the giant magnon can be thought of as a solution of the light-cone 
gauge-fixed string sigma-model described by a solitonic profile y = y{a — vt), where 
y is one of the fields parametrizing the five-sphere and v is the velocity of the soliton"^. 
In the infinite P+ limit this soliton exhibits the dispersion relation (2.10), where p 
coincides, in fact, with the total world-sheet momentum p^g carried by the soliton. 
Owing to the same form of the dispersion relation in the dual gauge theory, this 
gives a reason to call this soliton a "giant magnon" [37]. For our further discussion 
it is important to realize that if, instead of taking the field y from the five-sphere, we 
would make a solitonic ansatz z = z{(J — vt), where z is one of the fields parametrizing 
AdSs, we would find no solutions exhibiting the dispersion (2.10). As we will now 
show, in the mirror theory the situation is reversed: this time the giant magnon 
propagates in the AdS part, while there is no soliton solution associated to the five- 
sphere. 

Take the string Lagrangian (9.1) in the gauge a = 1 and put all the fields to 
zero except a single excitation z from AdSs. Upon making the double Wick rotation 
(2.16), the corresponding mirror action can be written as follows 



s = w''djdf -i + vl±£!^i:!±(l±i!)£!Urd5df^, (2.17) 



1 + 

Here r is an integration bound for a and z = dfZ , z' = d^z. Although our goal is to 
identify the mirror magnon configuration in the decompactification limit, i.e. when 
r — * oo, for the moment we prefer to keep r finite. 

To construct a one-soliton solution of the equations of motions corresponding to 
the action (2.17), we make the following ansatz 

z = z{a — vf) . 

Our further discussion follows closely [38]. Plugging the ansatz into (2.17), we obtain 
the reduced Lagrangian, Lred = Lred{z, z'), which describes a one-particle mechanical 
system with a treated as a time variable. Introducing the canonical momentum vr 
conjugate to z, we construct the corresponding reduced Hamiltonian 



i^rcd = 'JTz' - L 



red ' 



which is a conserved quantity with respect to time a. Fixing Hj-^d = 1 — u, where u 
is a constant, we get the following equation to determine the solitonic profile 



(^0^= /^f f, - (2.18) 



*See [38] and appendix 9.1 for more details. 
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The minimal value of z corresponds to the point where the derivative of z vanishes, 
while the maximum value is achieved at the point where the derivative diverges 



u;^ V f 



1 

1 , V <ijJ <\. 



The range of a is determined from the equation 

\z' 



da 







— = — f 2 E ^arcsin(zv^Ci;2/ (1 — uj'^)),!]^ 



where we have introduced r/ = ^jE^- Here E stands for the elliptic integral of the 
second kind. We see that the range of a tends to infinity when u ^ 1. Thus, u ^ 1 
corresponds to taking the decompactification limit. 

The density of the world-sheet Hamiltonian is given by 

where Pz = ^ is the momentum conjugate to z with respect to time f. For our 
solution in the limiting case = 1 we find 

v\z\ 

Pz -- 



^l-v\l + z^) ■ 
The energy of the soliton is then 



H = g da Ti = 2g I —— Ti = 2g arcsinh 



To find the dispersion relation, we also need to compute the world-sheet momen- 
tum p^s, the latter coincides with the momentum p of the mirror magnon considered 
as a point particle. It is given by 



P=Pws = - / dap^z' = 2 dz\pz\ = 2 — -- — . 

Finally, eliminating v from the expressions for H and p we find the following disper- 
sion relation 

\p\ 

H = 2q arcsinh — . 

To consider the semi-classical limit g ^ oo, one has to rescale the time as f —>■ f/g 
so that the energy H — > gH will be naturally measured in units of 1/g. Under this 
rescaling the momentum p scales as well, so that the dispersion relation takes the 
form 

\p\ 

i7 = 2 arcsinh ^ , (2.19) 

which is precisely the previously announced expression (2.14) for the energy of the 
mirror magnon. 
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3. Mirror S-matrix and supersymmetry algebra 

The S-matrix in field theory can be obtained from four-point correlation functions by 
using the LSZ reduction formula. Since the correlation functions can be computed 
by means of the Wick rotation, it is natural to expect that the mirror S-matrix is 
related to the original one by the same analytic continuation 

S{P1,P2) = S{pi,P2) , (3.1) 

where we replace Pi in the original S-matrix by Pi by using formulas (2.12). Just as 
the original S-matrix, the resulting mirror S-matrix should satisfy the Yang-Baxter 
equation, unitarity, physical unitarity, and crossing relations for real pk- 

On the other hand, the original S-matrix is su(2|2) ©5u(2|2) invariant and the 
states of the light-cone gauge-fixed AdSs x S^ string theory carry unitary represen- 
tations of the symmetry algebra su(2|2) © su(2|2). Therefore, if the relation (3.1) is 
correct then the mirror S-matrix should possess the same symmetry, and the states 
of the mirror theory also should carry unitary representations of su(2|2) ©su(2|2). 
This indicates that there should exist a way to implement the double Wick rotation 
on the symmetry algebra level, and that is what we discuss in this section. 

3.1 Double Wick rotation for fermions 

It is obvious that the double Wick rotation preserves the bosonic symmetry SU(2)^. 
To understand what happens with the supersymmetry generators it is instructive 
to apply the double Wick rotation to fermions. We consider the quadratic part of 
the light-cone gauge-fixed Green-Schwarz action depending on the fermions t] in the 
form given in [57] 

^ = WaVa - ^ {vaV5-a " vlv'l-a) " vlVa = WaVa " ^ . (3.2) 

Here, a = 1, 2, 3, 4, and we set k = 1 and rescale a in the action from [57] so that A 
disappears. 

Computing again the partition function of the model and using the path integral 
representation, we get 

Z{R,L) = j Vr]^Vr]e^o''d^fo 'i-(~vUa~n) _ (3 3) 

We note that fermionic variables here are anti-periodic in the time direction: 

rj{T + R) = -r]{T). 

Would fermions be periodic in the time direction, the corresponding path integral 
would coincide with Witten's index Tr(— l)^e~'^^, where F is the fermion number 
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[58] . Since in the mirror model r plays the role of the spatial direction, the mirror 
fermions are always anti-periodic in the spacial direction of the mirror model. On 
the other hand, the periodicity condition in the time direction of the mirror model 
coincides with a fermion periodicity condition in the spacial direction of the original 
model. In particular, if the fermions of the original model are periodic then the 
partition function of the original model is equal to the Witten's index of the mirror 
model. 

After the first Wick rotation the Lagrangian takes the form 

^ = -vlVa - ^ {vaVLa ' vlv'l-a) " vlVa ■ (3.4) 

Note that the fermions in this Euclidean action are not anymore hermitian conjugate 
to each other. 

Let us now perform the following change of the fermionic variables 

= ^ (i^Ua - V^a) , = ^ W + V^5-a) • (3.5) 

Computing (3.4), we get 

^ = -i^i^'a - I (^a^5-a - V^l^^a) " i^l^a • (3.6) 

It is the same action as (3.4) after the interchange t ^ a and ip ^ r], and this shows 
that the double Wick rotation should be accompanied by the change of variables (3.5). 
Note, that in terms of ip^s the supersymmetry algebra has the standard form with the 
usual unitarity condition. Thus, we expect that the supersymmetry generators will 
be linear combinations of the original ones. One may assume that in the interacting 
theory (beyond the quadratic level) one would take the same linear combinations. 

To summarize, the consideration above seems to indicate that the symmetry alge- 
bra of the mirror theory should correspond to a different real slice of the complexified 
su(2|2) ©su(2|2) algebra. Moreover, one might expect that the unitary representa- 
tion of the AdSs x string model could be chosen in such a way that its analytic 
continuation by means of formulae (2.12) would produce a unitary representation of 
the mirror model. 

3.2 Changing the basis of supersymmetry generators 

Let us recall that the centrally extended su(2|2) algebra consists of the bosonic 
rotation generators hj' , Rq,^, the supersymmetry generators Qo", Ql", and three 
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central elements H, C and C^. The algebra relations are 

[L,^ J^] = -6:3' + , [Rj, 3^] = -613^^ + ^6^J^ , 

{Q„^ Q/} = e.pe^" C , {Ql", qJ/^} = e.,e'^^ . (3.7) 

Here in the first two lines we indicate how the indices c and 7 of any Lie algebra 
generator transform under the action of L^^ and R^^- For the AdSs x string 
model the supersymmetry generators Qq,*^ and Q^", and the central elements C and 
are hermitian conjugate to each other: (Qa")^ = Ql"- The central element H is 
hermitian and is identified with the world-sheet light-cone Hamiltonian. It was shown 
in [10] that the central elements C and are expressed through the world-sheet 
momentum P as follows 

C = % (e'^ - l)e^^^ Ct = -%(e-^P-l)e-^^^ g = ^. (3.8) 

The phase ^ is an arbitrary function of the central elements, and reflects the ob- 
vious U(l) automorphism of the algebra (3.7): Q e*^Q , C e^*^C. In our 
previous paper [47] we fixed the phase ^ to be zero to match the gauge theory spin 
chain convention [9] and to simplify the comparison with the explicit string theory 
computation of the S-matrix performed in [26]. As we will see in a moment, if we 
want to implement the double Wick rotation under which P — > iH , H iP on 
the algebra level then we should choose ^ = — P/4. This choice makes the central 
elements C and to be hermitian and equal to each other^ 

C = Ct = -^sin-. (3.9) 
2 

As we discussed above, the symmetry algebra of the mirror theory should cor- 
respond to a different real slice of the complexified su(2|2) ©su(2|2) algebra. This 
means that we should give up the hermiticity condition for the algebra generators and 
consider a linear transformation of the generators which is an automorphism of the 
complexified 5u(2|2) ©5u(2|2) algebra. The transformation (3.5) suggests to consider 
the following change of the supersymmetry generators which manifestly preserves 
the bosonic SU(2)^ symmetry 

Q^'^=±_{Q^--^ Qt7 g^^) , Ql^ = ^ (QI" - z e"^ Q/ e,„) . (3.10) 



possibility of this choice was noticed in [10]. 
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Then, by using the commutation relations (3.7), we find 

{Q«^ Q;^} = 5,"R/ + 5f L," + '-515^^ (C + Ct) , (3.11) 

{Qa^ Q/} = e^pe'^' ^(C - Ct + zH) , {Qt", Qt/^} = e^.e'^^ 1{C^-C + lU) . 

Now we see that if we want to interpret the change of the supersymmetry generators 
as a resuh of the double Wick rotation then we should choose the central elements 
0,0^ to be of the form (3.9) because with this choice the algebra relations (3.11) 
take the form 

{Q«", Ql^} = S^^J + 5f L,-^ - 2zg sin ^ , 

{Qa^ Q/} = eo^^e^' -H , {Qt", qI^} = 6,,6"^ -H , 

and performing the analytic continuation 

we obtain the mirror algebra 

{Qa^ Ql^} = ^6 R/ + €W + 9 StS^ sinh ^ , 

^ ^ (3.13) 

Note that after the analytic continuation has been done we can impose on the new 
supersymmetry generators Q and new central elements H , P the same hermiticity 
condition as was assumed for the original generators. It is also clear that the algebra 
(3.13) implies the mirror dispersion relation (2.11). 

3.3 Mirror S-matrix 

The symmetric choice of the central charges (3.9) differs from the one we made in [47]. 
The S-matrix corresponding to the symmetric choice (3.9) coincides, however, with 
the string S-matrix in [47]. Indeed, this choice simply corresponds to multiplication 
of Q and Q} by e~*^/^ and e*^/^, respectively, which apparently does not change 
the invariance condition for the S-matrix. On the other hand, the string S-matrix 
also depends on the parameters ?7's which refiect the freedom in the choice of a 
basis of two-particle states. This freedom was partially fixed in [47] by requiring the 
string S-matrix to satisfy the standard Yang-Baxter equation. This still allows one 
to change the basis of one-particle states, or, in other words to change the basis of 
the fundamental representation of su(2|2). We will see that the requirement that the 
representation remains unitary after the analytic continuation fixes the parameters 
?7's basically uniquely. 
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To this end, we compute the action of the generators Q and Q''' on the funda- 
mental representation of su(2|2), see [9, 11, 47] for details. Starting with 



(3.14) 



we get 



(3.15) 



where 



a=—^{a + ic), b = —={b + id) , c=—^{c + ia), d = —^{d + ib) . (3.16) 
v2 v2 v2 v2 

and Qa", Qjj" are defined by eqs.(3.10). The unitarity of the representation after 
the analytic continuation requires 

(c + ia) = b* -id* . (3.17) 

The parameters of the original unitary representation before the analytic continua- 
tion are given by 



igx 




igx^ 






2 










igx^ 


g-«(?+'P) 




2 







1 jigx - igx+_^i(^-^) 
X' " 



d=d ^ e-^^--> , c = --^ 2 ' 



1 l igx ~igx^ 



(3.18) 



where ^ ~ p and (f ^ p are real, and the parameters x^ satisfy the following complex 
conjugation rule 

(x+)* = x~ . (3.19) 
After the analytic continuation, ^ , (f and p become imaginary (so that p is real) and 

(x+)* = — . (3.20) 
x~ 

Taking this into account and computing (3.17), we find that the analytically contin- 
ued representation is unitary for any choice of ^ if 




e^*^ = W— = e5P. (3.21) 



This means that the S-matrix which is unitary for real p and real p is obtained from 
the string S-matrix, see eq. (8.7) in [47], by choosing 

Vi = ?7(pi)e^^^ , r]2 = r]{p2) , r/i = r]{pi) , ^2 = r]{P2)e^^^ , (3.22) 
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where we have introduced 

v{p) = e^^^/ix~{p) — ix+{p). (3.23) 
Up to a scalar factor the S-matrix reads as [47] 



S{p, , P2) = "^(eI® El + El ® El + El ® El + El ® El 

xl - 1J1V2 ^ 



^ (eI (E) EI + EI(S) El - Ef (g) El ~ El ® Ef) 
mn-) \ 1 



(x^ -xWx^x^ ~ x^x^) mm 
(sf (g) El + Ej (g) Ej + El (E) Ej + E 

(.xr - xt){x^ - xt)ixT^+x+) ^ ^4 ^ ^4 ^ ^3 _ ^4 ^ ^3 _ ^3 ^ 

[X^ X2 )\Xi X2 X1X2) 
X2 -X^ TJl ^ ^3 _^ ^1 ^ ^4 _^ ^ ^3 _^ ^2 ^ ^ 

- x^ I'll V 

'e^I ® i^i^ + i;4 ® + £;f «) ^2 + ® i;. 



X+ - X^ 772 



+ I 



Xi — X2 ??2 
{x^ — xl\){x2 ^ X2){xl\ — X2) 

{^1 - - xr[x2)mm 



(e^ ®EI + EI®e1-E^2® El - El ® E^ 



- + +/ - _ — +W-. _ — _ I ^3 w i^4 JJ/4 w -t^3 - i^4 w J^3 - i^3 w ij| 

x-^ X2 \Xi X2 )(^ x^ ^2 ) 

Xt-Xiri2f^,^ £l ^ £4 ^ £l _^ £3 ^ ^2 _^ £4 ^ £2 



c+ - X2 m 



x\ f]2 



E\ ® El + E\ ®Ef + El®El + Ej® E^'^ , (3.24) 



where with i, j = 1, . . . ,4 are the standard 4x4 matrix unities, see appendix A 
of [47] for notations. 

With the choice (3.22) the S-matrix (3.24) satisfies the Yang-Baxter equation 
and it is unitary for real p's. The analytically continued S-matrix S{pi,p2) is then 
obtained from (3.24) by simply substituting 

1 



p— ^2iarcsinh — yl+p^, (3.25) 
c.f. section 2.2. One can verify that this matrix is also unitary for real p's: 

S{p,,p2)S^{pi,P2) = l. (3.26) 



The only subtlety here is that the string S-matrix also depends on a scalar factor, 
which has been omitted so far. Thus, one should separately check that this factor 
remains unitary after the analytic continuation. This will be discussed in section 5.2. 
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An exact relation between the S-matrix, 5*^^^, found in [47] and the S-matrix 
(3.24) is given by the following transformation® 



where G{p) = diag(l, 1, e*?, e**). It is amusing to note that a similar transformation 
has been recently introduced in [60], but with a very different motivation. Namely, 
as was shown in [60], the graded version of 5(^1,^2) coincides with the Shastry R- 
matrix [61] for the one- dimensional Hubbard model [62]- [64]. In section 4 we will 
give another interesting interpretation to our choice (3.22) which is based on the 
requirement of generalized unitarity. We will also show there that this choice of 
?7's makes the S-matrix (3.24) and, therefore, the Shastry R-matrix a meromorphic 
function on the ^-torus. 

To summarize, in order to have a unified description of the symmetry algebra 
of the AdSs x S^ light-cone gauge-fixed string theory and its mirror sigma-model 
we should make the symmetric choice of the central charges (3.9), and choose the 
fundamental representation of the centrally-extended su(2|2) with the parameters 
a, b, c, d given by 



Taking into account (3.16), (3.19) and (3.20), it is easy to check that both the original 
and the mirror (analytically-continued) representations are unitary with respect to 
their own reality conditions. Let us stress that the parameters a, b, c, d have the 
same dependence on in the original and mirror theories. We simply regard as 
functions of p in the original model, and as functions of p in the mirror one. 

3.4 Hopf algebra structure 

Formulas (3.27) define how the algebra generators of the original and mirror theories 
act on one-particle states of the theory. We also need to know their action on an 
arbitrary multi-particle state. The simplest way to have a unified description of their 
action is to use the Hopf algebra structure of the unitary graded associative algebra 
A generated by the even rotation generators 'Lj' , Ra^, the odd supersymmetry 
generators Qa", Q]j° and two central elements H and P subject to the algebra 
relations (3.7) with the central elements C and expressed through the world-sheet 
momentum P by the formula (3.9). We will be using the Hopf algebra introduced 
in [47] which is basically equivalent to the Hopf algebras discussed in [65], see also 

^The finite-size correction to the dispersion relation found in [32] involves the coeflicients oi, 02 
and ag of 5*^^^ (see [47] for notation) which are unaffected by this transformation. 



S{P1.P2) = Gi(pi)G2(p2)5^^"(Pl,P2)G'i(j9i)^^G2(p2) 




(3.27) 
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[66] for further discussion of algebraic properties of the centrally extended su(2|2) 
algebra. 

Let us recall that the unit, e : ^ ^ C, is defined by 

e(id) = l, e(J) = 0, e(Q) = 0, e(Qt) = , (3.28) 

and the co-product is given by the following formulas^ 

A(J) = J ® id + id ® J , 
A(Q„") = ® (Fl^ + e-'P/^ ® Q^-^ , (3.29) 

where J is any even generator. Here we use the graded tensor product, that is for 
any algebra elements a, 6, c, d 

(a ®h){c®d) = {-\)<^'^'^^\ac ® bd), 

where e(a) = if a is an even element, and e(a) = —1 if a is an odd element of the 
algebra A. 

It is interesting to note that the antipode S is trivial for any algebra element, 
that is 

3(3) = - J , SiQ) = -Q , 5(Qt) = -Qt . (3.30) 

This action of the antipode arises for the symmetric choice (3.9) of the central ele- 
ments C and only. 

The co-product is obviously compatible with the hermiticity conditions one im- 
poses on the algebra generators in the AdSs x string theory, and this ensures 
that the tensor product of two unitary representations is unitary. To check if the 
co-product is also compatible with the hermiticity conditions one imposes on the 
algebra generators of the mirror model we compute the co-product action on the 
supersymmetry generators Q , 

HQa") = Qa" ® cosh i^ — j + cosh (^ — j ® Q^," 

+ie^'^Qlfesa ® sinh j - i sinh j ® e^'^Ql^^ega , (3.31) 
A(Ql") = Ql" ® cosh (I) + cosh (I) ® Qt" 

-^e"'Q/eda ® sinh + ,sinh (^^j ^ e^'^Q/e^, . 

''To derive these expressions from the ones given in [47] one should rescale the supersymmetry 
generators in [47] by e**^/^. 
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Since in the mirror theory H is hermitian, the co-product is also compatible with the 
hermiticity conditions of the mirror theory. This guarantees that an su(2|2)-invariant 
S-matrix can be always chosen to be unitary. 

The co-product (3.31) can be used to find the commutation relations of the 
supersymmetry generators with the Zamolodchikov-Faddeev (ZF) operators A(j)) 
and A^{p) which create asymptotic states of the mirror model. The relations can 
be then used to determine the antiparticle representation, and to derive the crossing 
relation following the steps in [47]. A simple computation gives 

Q„Mt(p) = At(p)^„'^cosh + cosh (^)At(p)SQ,'^ (3.32) 
+ iA^ip){e''''^/esa) sinh (^) - tA^ip) sinh (^^^(e'^'^Qt^e^^) , 

Qi"At(p) = At(p)^,°cosh(|) +cosh(|)At(p)SQt" (3.33) 

- ^At(p)(e°^^/erf,) sinh +,sinh At(p)S(e"^Q5'^e,,) , 

where ^q," and are the matrices of the symmetry algebra structure constants 
corresponding to the fundamental representation (3.27) and S = diag(l, 1, —1, —1). 

As was already noted, the unitarity of the mirror S-matrix can be, however, 
broken by a scalar factor. In section 5 we show that the physical unitarity of the 
mirror S-matrix (the scalar factor) follows from the crossing relations. 

4. Double Wick rotation and the rapidity torus 
4.1 The rapidity torus 

The universal cover of the parameter space describing the fundamental representa- 
tion of the centrally extended su(2|2) algebra is an elliptic curve [12]. Indeed, the 
dispersion formula 

H^-4g^sm^^ = l, (4.1) 

which originates from the relation between the central charges of the fundamental 
representation, can be naturally uniformized in terms of Jacobi elliptic functions 

P 

p = 2a.mz, sin - = sn(z, /c) , H = dn{z,k), (4.2) 

where we introduced the elliptic modulus® k = —4:g'^ = — < 0. The corresponding 
elliptic curve (the torus) has two periods 2uji and 2u!2, the first one is real and the 

*Our convention for the elliptic modulus is the same as accepted in the Mathematica program, 
e.g., sn(z, fc) = JacobiSN[z, fc]. Since the modulus is kept the same throughout the paper we will 
often indicate only the z-dependence of Jacobi elliptic functions. 
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second one is imaginary 

2a;i = 4K{k) , 2uj2 = 4zK(l - k) - 4K{k) , 

where K(/c) stands for the complete elhptic integral of the first kind. The dispersion 
relation is obviously invariant under the shifts of z by 2uJi and 2uj2- The torus 
parametrized by the complex variable z is an analog of the rapidity plane in two- 
dimensional relativistic models. 

In this parametrization the real z-axis can be called the physical one for the 
original string theory, because for real values of z the energy is positive and the 
momentum is real due to 

1 <dn{z,k) < Vt , zeR, 

where k' = 1 — k is the complementary modulus. 

We further note that the representation parameters x^, which are subject to the 
following constraint 

,1 _ 1 2z 

X H — — — X = — , 

x+ x~ g 

are expressed in terms of Jacobi elliptic functions as 

= l/'S^i.Vl + dnz). (4.4) 
2g\s\iz J 

This form of x^ follows from the requirement that for real values of z the absolute 
values of x^ are greater than unity: > 1 if 2 G M. Note also that for real values 
of z we have Im(x"'") > and Im(x^) < . 

Since both the dispersion relation and the parameters x^ are periodic with the 
period c^i, the range of the variable Re z can be restricted to the interval from —uJi/2 
to uJi/2 which corresponds to — tt < p < tt. 

Postponing an extensive discussion of the bound states till section 7, we note here 
that the latter problem requires consideration of complex values of particle momenta. 
According to eq.(4.2), a rectangle —uJi/2 < Yie^z) < uJi/2 ; —uj2/2i < \m.{z) < uJ2/2i 
is mapped one-to-one onto the complex p-plane. By this reason, it is tempting to 
call this rectangle by the physical region in the complex z-plane,^ and, therefore, to 
restrict the allowed values of the 2;-coordinates of the particles forming a bound state 
by this region. An advantage of adopting such a choice is that all the bound states 

^In relativistic field theories treated in terms of the rapidity ^ 62 — Oi, the physical region is 
defined as a strip < Im0 < tt and it incorporates the bound states. Correspondingly, the physical 
region of an individual particle is ImO G (— 7r/2,7r/2) and it covers the complex p-plane (with a 
cut) through the relation p = sinh 9. 



(4.3) 
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Figure 1: On the left figure the torus is divided by the curves = 1 and = 1 into 
four non-intersecting regions. The middle figure represents the torus divided by the curves 
Im(x^) = 1 and Im(a;~) = 1, also in four regions. The right figure contains all the curves 
of interest. 

would have positive energy. We will see, however, that this is not the only option, 
and there are other two regions in the complex 2;-plane which could equally deserve 
the name ^^physicaT . As it will become clear later on, counting the degeneracy of 
the bound states drastically depends on the choice of a physical region. 

Each solution of eq.(4.3) corresponds to a point of the half-torus, i.e. of the 
rectangle^" —LOxjl < Re(z) < uJi/2; — 3(^2/4^ < Im(z) < 5u;2/4z. In what follows we 
will be loosely referring to this rectangle as the torus. The torus covers the complex 
p-plane twice. Since the space of solutions of eq.(4.3) is mapped one-to-one on the 
torus, the latter could be also chosen as the physical region. Such a choice is however 
problematic because half of all the states would have negative energy, i.e. the region 
would contain both particles and anti-particles, as well as bound states and anti- 
bound states. We point out, however, that there exist positive energy solutions of 
the bound state equations with some of the particles falling outside of the rectangle 
—uJi/2 < Re{z) < uJi/2; —uJ2/2i < lm{z) < uJ2/2i that covers the complex p-plane 
once. 

Constraint (4.3) implies that if a pair {x^,x~) satisfies it then x~), 
[x^ , l/x~) and also do. Each of these four pairs corresponds to a 

different point on the torus. Taking into account that for any complex number w 
if \w\ > 1 then < 1, and if Im(w;) > then lm{l/w) < 0, one can divide the 

torus into four non-intersecting regions in the following two natural ways, see Fig.l: 

^"We made slightly asymmetric choice for Im(z) to achieve better visual clarity. 
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g=1/2 g=1 g=50 

Figure 2: Divisions of the torus by the curves |x^| = 1 (upper figures) and by the curves 
Imx^ = (lower figures) for g = 1/2, g = 1 and g = 50. The red curves are = 1, and 
the pink ones are \x~^\ = 1. The coordinates x and y are the rescaled real and imaginary 
parts of z: x = Re(^z), y = Re(^z). In the limit 5 — > oo the curves = 1 and 
Irax^ = are related by the shift z ^ z + 

• {\x^\ > 1}, < 1}, < > 1} and > < 1}; the 
division is done by the curves = 1. 

• {Im(a;=^) > 0}, {Im(a;=^) < 0}, {Im(a;+) > 0, Im(x") < 0} and {Im(a;+) < 
, Im(a;^) > 0}; the division is done by the curves Im(a;^) = 0. 

The shape of the regions depends on the value of the coupling constant g, see 
Fig. 2. Quite remarkably, in the strong coupling limit g ^ 00 two divisions of the 



-24- 



torus produced by the red (|a;^| = 1) and green (Im(a;^) = 0) curves become related 
to each other through a global shift by C1J2/2. 

There are eight special points on the torus where the curves \x^\ = 1 intersect 
with the curves Im(x^) = 0, see Fig.l. These points are z = i^cji + ^^^^ci;2 , n = 
—2, —1, 0, 1. It is known [18] that these points are the branch points of the one-loop 
correction [16] to the dressing phase. It is unclear, however, if they remain the branch 
points of the exact dressing phase. One could try to use the integral representation 
[52] of the BES dressing phase [14] to understand this issue. In fact, all currently 
available representations for the dressing phase are defined for > 1, and this is 
another reason to figure out the location of the curves = 1 on the 2-torus. 

Both divisions play an important role in the analysis of the bound states of 
string and mirror theories. To understand the meaning of the equations = 1 
and Im(x^) = 0, it is convenient to use another parameter u which is similar to 
the rapidity parameter of the Heisenberg spin chain. In terms of it is defined as 
follows 

M = x+ + — = x^H h-. (4.5) 

g X g 

By using eqs.(4.5) and (4.4), one can express the rapidity m as a meromorphic function 
on the torus 

ciazdnz , , . 

u = . 4.6 

g snz 

It is not difficult to check that the eight special points on the torus are mapped 
onto the four points on the M-plane with coordinates u = ±2 ± -, while the points 
z = ±uJi/2 are mapped to m = 0, and the points z = ±uJi/2 + UJ2/2 ±iO are mapped 

to M = ±00 ± too. 

A special role of the points u = ±2 ± - can be also understood by expressing x^ 
in terms of u 



X 



X 




(4.7) 



Thus, on the u-plane there are four branch points with coordinates u = ±2 ± - 
corresponding to x^ = ±1 and Im(x^) = 0. Therefore, we can naturally choose the 
cuts either connecting the points — 2 ± ^ and 2 ± ^ , or going from ±00 to ±2 ± ^ 
along the horizontal lines. Let us determine what values of x^ correspond to the 
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Figure 3: On the left figure the upper and lower curves correspond to = 1 + 
and \x'\ = 1 + 0, respectively. The map z u{z) folds each of these curves onto the 
corresponding cut on the u-plane. 



lines u 



Mo ± ^ with uq real. We see that 



Uo = X^ + —, 



Uq = X + 



X~ 



X 



X 



^ ( Mo ± \/mo - 4 ) , if M = Mo - - , 
^ j^Mo ± y^Mg - 4j , ifM = Mo + ^. 



It is clear that points x^ and 1/x^ of the complex x^-plane correspond to the same 
point M of the M-plane. Then, the points of the circle \x~^\ = 1 map to points m in the 
interval [—2 — - , 2 — -] , while the points of \x~ \ = 1 correspond to m G [—2 + - , 2 + -]. 
On the other hand, the points of the lines Im(x"'") = and Im(x~) = correspond to 
points M outside the intervals [—2 — - , 2 — -] and [—2 + - , 2 + -], respectively. Note 
also that if one chooses a definite sign in eq.(4.7) then the interval [—2 + - , 2 + -] 
maps onto a half of a unit circle in the x^-plane. One has to use both signs to cover 
the unit circles = 1 and real lines Im(a;^) = 0. 

To determine the location of the upper and lower edges of the M-plane cuts 
[— 2 + ^ , 2 + ^] on the x^-planes, we introduce a small real parameter e and write 



x^ = e'e'^P 



„±i 



1x^1=6% Im(a; ) = e''siny9, m~2cosv9+ — |-2iesin(y9. (4.8) 



We see that the upper edges [—2 + ^ + zO , 2 + ^ + zO] are mapped either outside the 
upper halves or inside the lower halves of the circles = 1, and the lower edges 
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[—2 =p - — iO , 2 =p - — iO] are mapped either outside the lower halves or inside the 
upper halves of the circles = 1, and vice verse: 



- + zO,2=F- + zO <^ ' ) ,[ 

g ' ^ \ x± = 1 -0, Im < 



t ^ « f a;^ = 1 + 0, Im(x=^) < 

g ' g ^ \ = 1 -0, Im(x±) > 



As we discussed above, the 2;-torus can be divided into four non-intersecting re- 
gions by the curves | = 1. Now it is easy to show that each of the regions is mapped 
one-to-one onto the w-plane with the two cuts. Let us consider for definiteness the 
region with > 1. Then, according to the discussion above, the boundaries of the 
region with \x'^\ = 1 + 0, Im(x^) > and {x'^l = 1 + 0, Im(x'*") < are mapped 
onto the upper and lower edges of the cut [—2 — - , 2 — -] in the w-plane, respectively. 
In the same way the boundary of the region with |a;~ | = 1 is mapped onto the upper 
and lower edges of the cut [—2 + ^ , 2 + ^], see Fig.3. 

Another way to understand how different copies of the M-plane are glued together 
is to consider any of the curves = 1 and shift its variable z by a small positive 

e in the imaginary direction. For the image of the corresponding shifted curve on the 
M-plane one obtains 

1 /du\ 

Im m(z + ze) = +- + e Re — + . . . , (4.11) 



g \dz 

where Re (|^) is computed along |a:^| = 1. Further analysis shows that along any 
of the curves = 1 the expression Re (||) is positive for — ^ < Re 2; < ^ and 
negative otherwise. This determines how the edges of the cuts = 1 are mapped 
onto the edges of the corresponding cuts on the w-plane (see Fig.3 for an example ). 

To summarize, any region confined between the curves \x^ \ = 1 is mapped under 
z u{z) onto a single copy of the u-plane with a point at infinity added, i.e. onto the 
Riemann sphere. Extended to the whole torus, this map defines a four-fold covering 
of the Riemann sphere by the torus which has eight ramification points-}^ a generic 
point on the M-plane has four images belonging to the four regions. There are two 
cuts on each copy of the w-plane 

1) [-2 + i/g,2 + i/g] 

2) [-2-i/g,2-i/g] 

which are images of the curves = 1 and = 1, respectively. 

In the same way we can determine the images of the upper and lower edges of 
the M-plane cuts (—00, — 2 + ^] , [2 + ^, 00) on the x^-planes. We again introduce a 



11 



In agreement with the Riemann- Hurwitz formula. 
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Figure 4: Four copies of the u-plane (the Riemann sphere) glued together through the 
cuts to produce the torus of the kinematical variable z. We indicated four branch points 
Bi^2 and C1.2 which are images of those on Fig. 3. 

small real parameter e and write 

\ i 1 

= re**^ , |a;^| = \r\ , Im(x^) ^re, u ^ r -\ — =]= — h ie(r ) . (4.12) 

r g r 

We see that the upper edges (—00, — 2 =1= ^ + iO] , [2 =1= ^ + iO, 00) are mapped either 
onto the upper edge of the intervals (—00, —1] , [1, 00) or the lower edge of the interval 
[— 1, 1], and the lower edges (—00, — 2 =1= ^ — zO] , [2 =f ^ — iO, 00) are mapped either 
onto the lower edge of the intervals (—00, —1] , [1, 00) or the upper edge of the interval 
[— 1, 1] of the real lines Im(x^) = 0, and vice verse: 

(-00, -2 =F - + iO] U [2 =F - + iO, 00) 
9 9 

(-00, -2 T - - iOl U [2 ^ - - iO, 00) 
9 9 

Again, dividing the 2;-torus into four non-intersecting regions by the curves Im(a;^) = 
0, we see that each of the regions also maps one-to-one onto the w-plane with the 
two cuts. This gives a different (but equivalent) four-fold covering of the Riemann 
sphere by the torus. 

When a point on the 2;-plane runs along the curve jx"*"! = 1 or = 1 its image 
covers the corresponding interval on the M-plane twice. To appreciate this fact, let 



Im(x^) 
Im(x^) 

Im(x^) 
Im(x^) 



< 1 

+0, < 1 
.±1 >1 



-0 



-0, 



(4.13) 
(4.14) 
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us note that if z is, e.g., on the curve \x^\ = 1 then the points 

., = -.±1 + 1 (4.15) 

are also on this curve. Indeed, since = x~^{z)x~ (z*), we have 

where we have used the properties of Jacobi elhptic functions under the shifts by 
quarter-periods. In the same way one finds that if z lies on a curve = 1 then 
the points z± belong to another copy of = 1 which is obtained from the original 
one by the shift by uj2- Finally, using the properties of the Jacobi elliptic functions 
it is easy to show that u{z±) = u{z), i.e. the points z and z± have one and the same 
image on the w-plane. 

It is clear that the half of the torus and, therefore, the complex p-plane is mapped 
onto the M-plane twice. The coordinate u is real for real z, and in this case we can 
easily express it in terms of p [5] 



= icot^Wl + Vgin2^. (4.16) 

In the limit g ^ the relation (4.16) turns to the one between the rapidity and 
momentum variables of the Heisenberg spin chain; the latter describes the gauge 
theory at the one-loop level. This supports an idea that the physical region could be 
identified with a single copy of the u-plane, namely the one which maps to the region 
> 1 of the 2;-torus. There are certain advantages of such a choice which we will 
discuss later on. The main disadvantage is, however, that the region \x^\ > 1 is not 
big enough to cover the whole complex p-plane. 

It is interesting to see what happens with our three candidates for the physical 
region in the limits g ^ oo and g ^ 0. In the limit g ^ oo the periods of the torus 
have the following behavior 

uji-^ , u;2 ^ — if 5-^00. (4.17) 

9 2^ 

To keep the range of lm{z) finite, we rescale z as z ^ z/{2g), and the momentum as 
P ~^ pI 9- Then the dispersion relation (4.1) takes the relativistic form — = 1, 
the variable z plays the role of 6 because p = sinh z, and we have 

• The torus degenerates to the strip with — tt < lin{z) < vr and —oo < Re(z) < oo 

• The half-torus corresponding to the complex p-plane degenerates to the strip 
with — 7r/2 < Im(2;) < 7r/2 and — oo < Iie{z) < oo 
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• The region |a; | > 1 corresponding to the complex w-plane degenerates to the 
strip with — 7r/2 < lm{z) < 7r/2 and —oo < Re(z) < oo 

We see that both the half-torus and the region > 1 degenerate to the physical 
strip of a relativistic field theory. 

In the limit g ^ the periods of the torus have the following behavior 

cui — > TT , uj2^2i\ogg if g . (4-18) 

We see that all the three regions degenerate into the strip with — 7r/2 < Re(z) < 7i/2 
and —oo < lm{z) < oo. The properties of the S-matrix arising in the limit g ^ 
will be discussed in appendix 9.2. 

4.2 Double Wick rotation 

The 2;-torus can be also used to describe the mirror model. Since we know the 
relation between p = 2 am z and the mirror momentum p, we can express p in terms 
of z. Indeed, the equality 

2 am z = 2i arcsinh— ^/l +jp (4.19) 

implies 

p= -idnz. (4.20) 
The energy in the mirror theory takes the form 

~ \/k' 

= 2arccoth^^ — . (4.21) 
anz 

The formulae above show that real values of z correspond to imaginary p. Now 
we would like to understand for which values of z the corresponding values of p are 
real. One can see that if we shift the variable z by UJ2/2, z ^ z + UJ2/2, that is if we 
write 

p=-tdn(z+'^,k)=Vk^—, (4.22) 
\ 2 / cn^; 

then for real values of the shifted variable z the corresponding values of p are real as 
well. We also recognize here a close analogy with the relativistic case - making the 
double Wick rotation corresponds to the shift by a quarter-period on the rapidity 
plane. The function cn(z, k) has zeroes at z = —^uJi and z = (and dn(2;, k) has 
poles dX z = + UJ2) and z = \{uJi + UJ2)) which explains the apparent absence 

of the periodicity in p. Thus, when the shifted variable z runs from — ^cji to ^c^i the 
momentum p monotonically increases from —00 to +00 and it passes though zero for 
^ = 0. 
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One further finds that the parameters are expressed in terms of the shifted 
parameter z of the mirror model as follows 



X = -' ' " ' 



fl + ,v^^y (4.23) 



\/—kdnz 

We can now find how x^ are expressed in terms of the mirror momentum. Indeed, 
since 

' cn z \ 2 k 



snzJ 1 — dn^ z ' 



we deduce from eq.(4.23) that 




x^ = 

This, of course, agrees with the formula (2.15). 

The variables x^ of the mirror theory obey a relation 

+ _ P-i 
X X = — 

p + i 

which implies that |a;"'"x~ | = 1 for p real. 

It is also not difficult to show that the dispersion relation in the mirror theory 
takes the form (2.11) 



_ k 1 

H = 2 arccoth- = 2 arccoth^ / 1 — = 2 arcsinh ^Jl +jp . 

This completes the proof that the double- Wick rotation corresponds to a shift of the 
z variable by a quarter of the imaginary period of the torus, and the real axes of the 
shifted z corresponds to real values of the momentum of the mirror theory. 

Finally, it is useful to express the rapidity u in terms of the shifted parameter z 
of the mirror model and p. We have 



u 



2cn{z + f,k) dn{z + f,k) _ 2iVk' dn {z + f , k) 
y/^sn{z + ^,k) y/^dn {z,k) 



Then one can check that the points z = ±uJi/2 ± iO are mapped to n = ±oo ± i oo. 
The coordinate u is real for real z, and in this case we can express it in terms of p 




]P g\i 1 + 



After having performed the shift, one can do various physically equivalent transformations of the 
shifted z-variable preserving the axes of real z. Particular useful examples of these transformations 

are z ^ z + ^ , z ^ -z + ^ , z ^ -2± ^. 
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Again, there are three choices of the physical region. It is the half-torus corre- 
sponding to the complex p-plane, the whole torus, and the region Im(a;^) < which 
is mapped onto the w-plane. The third choice is different from the one made for 
the string theory, and is motivated by the analysis of the bound states of the mirror 
model. 



5. S-matrix on elliptic curve 

5.1 Elliptic S-matrix and its properties 

The dispersion relation (4.1) is naturally parametrized by the elliptic curve. Without 
imposing the unitarity condition for the S-matrix, the phase rj in (3.23) can be chosen 
in an arbitrary way, for instance, T]{p) = 1. In the latter case, the S-matrix (3.24) is 
well defined on the elliptic curve but it is non-unitary. It is therefore tempting to as- 
sume that the unitary S-matrix also admits an analytic continuation into the complex 
2;-plane. To find such a continuation one has to resolve the branch cut ambiguities 
arising due to the ?7-factor in the S-matrix (3.24): r](p) = ei"^ ^yix~{p) — ix+{p). 

This can be done in the following way. First, we recall the elliptic parametrization 
(4.4) which gives 

vip) = \/ix~{p) — ix^{p) = ^—et*^™^ V 1 + dn 2; = 

= a/(1 + dn2;)(cn2; -|- isnz) . (5.1) 

V9 



Second, by using the following formulae (recall k = —Ag 



2^ 



2dn^| (cn f + i sn fdn f )^ 

1 + dnz = ; — ^ , cn z + i sn z— ^ ^ ^ 



l-fcsn^l' 1-A;sn4| 

relating elliptic functions to those of the half argument, we can resolve the branch 
cut ambiguities by means of the relation 

'r, I V7-^ \/2 dn f (cn f isn fdn f) , ^ , ^ 

e^^ ^ix~{p)-ix^{j>) = ^ i+%2s^4| = viz) (5.2) 

valid in the region — ^ < Re^ < ^ and —uj2/i < Imz < 002/1. Further, we notice 
that the non-local dependence of r/'s on the momentum of another particle enters as 
gfp _ gjamz Qj^f^^ therefore, can be naturally treated as e^^ = cnz + isr^z. 

Thus, we define an analytic continuation of the S-matrix onto the rapidity torus 
for each of the complex variables Zi and Z2 by means of eq.(3.24), where the variables 
?7i^2 and are given by 

7]i = 7]{zi){cnz2 + isnz2), V2 = viz2) , 

f/2 = 77(^2) (cn^i +isn 2:1) , 771 = 77(2:1). 
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In this way we completely resolved the branch cut ambiguities of the S- matrix (3.24) 
and defined it as a meromorphic function on the elliptic curve (for each 2;- variable) . It 
is remarkable to observe that such a continuation becomes possible due to additional 
phase factors, e^^, introduced in the previous section to guarantee unitarity of the 
mirror theory. 

Let us now analyze the basic properties of the elliptic S-matrix. One can check 
that it satisfies the Yang-Baxter equation and the usual unitarity requirement 

Si2{Zl,Z2)S2l{z2,Zi) =1. (5.3) 

Further, it obeys the generalized unitarity condition: 

^12(^!,Z2*)['5i2(Zl,^2)]^ = I. (5.4) 

Here "f" means hermitian conjugation. For zi and Z2 real the last condition reduces 
to the requirement of physical unitarity. In fact, one can see that the elliptic S-matrix 
is compatible with the generalized unitarity condition only due to our specific choice 
for the phase factors discussed above. Then, unitarity and generalized unitarity 
imply hermitian analyticity: S2i{z2, z^) = [Si2{zi, 22)]^ . 

Let us now compute monodromies of the S-matrix (3.24) over the real and imag- 
inary periods. We find 

S{zi + 2ui, Z2) = Si S{zi, Z2)T.i = T.2S{zi, ^2)^2 , , ^, 

(5.5) 

S{Zi + 2U2, Z2) = Si S{Zi, Z2)T.i = T.2S{Zi, Z2)'T^2 ■ 

Hence, the S-matrix exhibits the same monodromies over real and imaginary cycles 
and it is a periodic function on a double torus with periods and 4u;2- Here 
Si = S ® I and S2 = I ® S, where S is defined in section 3.4, and the S-matrix 
commutes with the product S ® S. Note that S is in the center of the group 
SU(2) X SU(2). 

Second, we establish the monodromy properties w.r.t. shifts by half-periods. 
Under the shift by the real half-period we get 

^(^1 + cc^i, Z2) ={y®Y) S{z^, Z2) {V-^ ® I) , (5.6) 

where V = diag (e i ,e ,e ^ ,e i j. 

The shift by the imaginary half-period corresponds to the crossing symmetry 
transformation [12]. To discuss it, we multiply the S-matrix (3.24) with a scalar 
factor So to produce the string S-matrix obeying crossing symmetry 

'5(^1,2:2) = So{Zi,Z2)S{zi,Z2) . (5.7) 

We then find that with a proper choice for 5*0(2:1, Z2) the string S-matrix exhibits the 
following crossing symmetry relations 

"^.-'Silizi, ^2)^1^12(^1 + 102, Z2)=I, WU^l, Z2)'^i-'Si2{Zl - CO2, Z2)=I, (5.8) 



33- 



and also 

^2)^l5i2(Zl, Z2 - ^2) = I , "^iSU^U Z2)'^r'5i2(^l, Z2+UJ2)=l. (5.9) 

Here ti denotes transposition in the first matrix space and is a constant^^ charge 
conjugation matrix 

where (T2 is the Pauh matrix. The compatibihty of eqs.(5.8) and (5.9) with (5.5) is 
guaranteed by the identity = "^"^ which is equivalent to = S. 

The crossing symmetry relations lead to the following equations for the scalar 
factor 5*0 [12] 

Sq{zi, Z2) Sq{zi + 102, Z2) = f{zi, Z2) , 5*0(21, Z2) Sq{zi, Z2 - 102) = f{zi, Z2) , (5.11) 
where the function / is expressed through as follows 



f{Zl,Z2) = 




(5.12) 



One can easily check that the function f{zi, Z2) obeys the following properties 

f{z2,Zi)f{Zi+U2,Z2) = 1 = f{z2,Zi)f{Zi,Z2+U2), f {Zi + U2, Z2 + UJ2) = f{Zi,Z2), 

which are, however, incompatible with the assumption that the scalar factor is an 
analytical function of Zi,Z2- 

Another important property of the string S-matrix (5.7) is that it remains in- 
variant under the simultaneous shift of zi and Z2 by 002'- 

S{zi + UJ2, Z2 + UJ2) = S{zi, Z2) . (5.13) 

This follows from the fact that both the S-matrix (3.24) and the scalar factor Sq are 
invariant under the shift. This property together with the crossing relations (5.8), 
(5.9) implies 

5*^'*^(Zi, ^2) = '^l'^2'5(2l, ^2)"^! ^'^2 ^ = '^l ^^2 ^'5(^1; ^2)*^! ^2 ; 

where ti and ^2 niean the transposition in the first and in the second matrix spaces, 
respectively. 

^■^This is in opposite to [47], where the charge conjugation matrix was found to depend on the 
sign of the particle momentum. This dependence is, in fact, spurious and it gets removed by a 
proper resohition of the branch cut ambiguities we propose here. 
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Assuming that the above-mentioned properties of the S-matrix (3.24) are shared 
by iS, we can now see that the string S-matrix allows one to define consistently an 
elliptic analog of the ZF algebra, i.e. 



A\iz^)Aliz2) = Al{z2)Al{z,)S,2{zi, z^) , 

where the creation and annihilation ZF operators are now functions of the complex 
variable z. In addition, away from the line zi = z^ we can impose the following 
relation between the creation and annihilation operators 

A^{zx)A\{z2) = Al{z2)S2i{z2, z^)Ai{zi) . (5.15) 

As usual, the absence of cubic and higher relations for the ZF operators is guaranteed 
by the Yang-Baxter equation for S. Furthermore, the validity of relations (5.14), 
(5.15) for all values of zi and Z2 is due to unitarity condition (5.3). 

Transposing the second equation in (5.14) in the first matrix space we get 

{A{{z,)Y^Al{z2) = Al{z2)Si}2{zuZ2){A\{zr)y^ , 

On the other hand, shifting in eq.(5.15) the variable zi by the imaginary half-period 
we obtain 

Al{Zi + LU2)AI{z2) = Al{z2)Si2{Zl + iU2, Z2y^ Ai{zi + UJ2) . 

Since the string S-matrix satisfies the crossing relation we see that the algebra struc- 
ture is compatible with the following identification 

A{z + U2) = "^-^ A\z)\ A\z-U2) = -A{zy'tf. (5.16) 

Analogously, we establish 

A{z-uj2) ='^A\zy , A\z + uj2) = -AizY'^-K (5.17) 

These relations together with the monodromy properties (5.5) of the S-matrix further 
imply 

A{z + 2ui) = J:A{z) , A\z + 2cji) = A\z)J: , 
A{z + 2002) = ^A{z) , A^{z + 2uj2) = A^{z)E . 

This means that the bosonic operators are unchanged under the shift around the 
torus while fermionic ones acquire the minus sign. Thus, the monodromy properties 
of the S-matrix imply the spin structure (— , — ) for the fermionic ZF operators. 
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Finally, the generalized unitarity condition (5.4) allows one to impose the fol- 
lowing hermiticity conditions on the ZF operators: 

[A'{z)]^ = AI{z*) for 0<|lm^|<^; 

n (5-18) 
[A\z)]^ = -4(z*) for ^ <\lmz\<^. 

The hermiticity condition for the ZF creation and annihilation operators in the anti- 
particle region uo2/2i < \lmz\ < 002/ i is compatible with the hermiticity condition 
for the ZF operators in the particle region < |Im2;| < 002/2% and the identifications 
(5.16) and (5.17). 

5.2 Unitarity of the scalar factor in mirror theory 

It is clear from the discussion above that the S-matrix of the mirror theory is obtained 
from the string S-matrix just by the shift of the z- variables by UJ2/2 

S{z,,Z2)=S{z, + ^,Z2 + ^). (5.19) 

The momentum of the mirror theory is expressed in terms of the variable z by 
eq.(4.20) and is real for real values of z, and the generalized unitarity of the mirror 
S-matrix in terms of the shifted coordinates z takes the usual form 

Siz^,Z2)]^ Sizlz*2)=I. (5.20) 

This just follows from the generalized unitarity of the string S-matrix and relation 
(5.13) which is a consequence of the crossing equations 



■>12[Zl,Z2 



- ^2l[Z2 ~ ~ ~2~ ~ '^2112:2 + —,Zj^ + —) — C>2l[Z2, 2i j • 



In fact, since both the S-matrix (3.24) and the scalar factor Sq satisfy the generalized 
unitarity condition and relation (5.13), the same holds for the mirror theory. 

It is of interest to understand how the dressing factor of the mirror theory trans- 
forms under the complex conjugation. To this end we recall that in the a = 
light-cone gauge^^ the scalar factor of the string S-matrix can be written in the form 
[18] 

+ 1 - +^ _ 

So{zi,Z2Y = s{zi,Z2)a{zi,Z2) , s{zi,Z2) = ^ (5.21) 

Here the gauge-independent dressing factor (t(zi, Z2) has the following structure [6] 

^ 00 00 

-\na{zi,Z2) = 6{zi,Z2) = ^ Q.s(fi') qr{zi)qs{z2) - qr{z2)qs{zi) , (5.22) 

r=2 s=r+l 

^'^It is easy to check that the additional a-dependent factor does not break any of the properties 
of the S-matrix. 
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where qr{z) = — (x^)^^*^'] are the local conserved charges. At any order of 

the perturbative expansion in powers ofl/g the sums in r and s define the convergent 
series for > 1 and \xf\ > 1. Thus, the S-matrix is by construction well-defined 
only in the region where > 1 and it should be analytically continued for other 
values of . 

The string theory dressing factor satisfies the generalized unitarity condition that 
follows from the fact that under the complex conjugation the variables x^ transform 
as = x'^{z*). In the mirror theory the variables x^ depend on the shifted 

coordinate z and, as a result, satisfy the following complex conjugation rule 



X [z ^ 

, 2 ' 



x^iz* + ^' 



By using this rule one can easily check that the factor 5(^1,^2) in (5.21) transforms 
under the complex conjugation as follows 



[s{zlzl)ts{z,,z,)=[^-^]\ (5.23) 

where xf = x^{zi + ^). Taking into account that the scalar factor of the mirror 
theory satisfies the generalized unitarity condition, we find the complex conjugation 
rule for the dressing factor of the mirror theory 



[aizl,z;)faiz^,Z2)=('^] . (5.24) 

In particular, for real values of z^s corresponding to real p^s the dressing factor of the 
mirror theory is not unitary. 

It is interesting to note that the scalar factor can be split into a product of two 
factors satisfying the generalized unitarity condition in both string and mirror theory 

^0(^1,^2)^ = - + X ^^(^1,^2). (5.25) 

Another interesting splitting is given by 




^0(^1,^2)' = X ( ^ I a{z,,Z2). (5.26) 

-U2 + J 

This splitting is useful for analyzing the bound state spectrum of the mirror model. 

Knowing the series representation for the dressing phase in the original the- 
ory [13], it is interesting to understand what is precisely the source of its unitarity 
breakdown in the mirror theory. 
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To clarify this issue, we recall that the dressing phase can be conveniently written 
in terms of a single function x{^iy^2) [18] 

- xixt, xt) + xix2,xi) + x{xt, X];) - xix2,x^) , 
which admits the following strong coupling expansion 



oo 

x{xi, X2) = X^'^Kxi, X2) ( 



n=0 

Here 



(0)/ N 1 X1X2-I, X1X2-I 

T '{Xi,X2) = log 

X2 X2 X1X2 

is the leading AFS factor [6]. The next-to-leading contribution is [16]: 



fl). X ly. y/Xl-l/y/x^ 1^. y/x^+l/y/x^ 
X {Xl,X2) = — ijl2 ijl2 

2n y/xi - y/x^ 2n .Jx{ + 



^Li, yxI+Vv^ + iLi^ v^-Vv^ . (5.27) 
27r ^Jx{ - ^Jx2 2% ^Jx{ + .Jx2 



All higher terms are rational functions of X\^X2 [13]. As we will now show, the 
unitarity breakdown of the dressing phase is due to the leading AFS contribution 
only, the Hernandez-Lopez term (5.27), as well as all higher order terms do not 
influence the unitarity condition. 

To simplify the notations in what follows we only consider the case of real 2;'s in 
the mirror theory. It is easy to see that the complex conjugate of the function 
is given by 



\)6^^{xt.x^')\ = -x^°\x^,x^) - ^^J^ ^ + {in - l)x^ - (^-^ - xj) logx^x 



2 • 



Using this formula for computing the leading value ^afs? we find that the contribution 
of non-logarithmic terms cancels out and we get 

^AFs = ^AFs + 9 ^ „ „^ (1 - x^x^) log 0:^X2 + g \ J {1 - xfx^) log X^X^ 

12 i-o ^ iX> 2 

- g ^ _ (1 - x^xt) logx^x^ + g \ _^ (1 - x^"X2 ) logx^^Xa" . 

Using identity (4.3), it is easy to show that all logarithmic terms are neatly combined 
to produce the following answer 

X1X2 ^ ^ 



'^1*^2 



eiFS = ^AFS + ^log^ , (5.28) 
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which coincides with the logarithmic form of eq.(5.24). 

Since we have shown that the shift of the phase under the complex conjugation 
occurs due to the leading contribution, all the higher order terms in the expansion of 
9 must be real functions. To convince oneself that this is indeed the case, we consider 
the next-to-leading term in the strong coupling expansion of 9. As was shown in [13], 
this term admits the following representation 

^HL = ^{qt - qt) - Hqt - q2) - i^{qi - qt) + iiqi - q2) ■ (5.29) 

Here the function il){q) is 

^lJ{q) = ^Li2(l - e^^) - ^Li^il - e'^+^-) - log(l - 6^"+^) + ^ , (5.30) 
where the variables are related to through 

e'»* = 2^. (5.31) 
Taking into account the conjugation rule in the mirror theory, eq.(3.20), we obtain 

(g^r = -g^-^. (5.32) 

Since 6'hl depends on the difference of two q's, the shift by tt arising upon the complex 
conjugation will cancel out. Thus, taking the complex conjugate we obtain 

^HL = i^iql - q2) - i^{qi - qt) - ^(qt - 9^) + i^iqt - qt) > (5.33) 

where the function ip{q) is defined as 

^(g) = ^U,{1 - e'^) - ;lLi2(l - e'^--) - log(l - e^"--) + ^ . (5.34) 

Taking into account the following transformation property of the dilogarithm func- 
tion 



Li2(l - e*^-^") = Li2(l - e^5+^"-2") = Li2(l - e^'?+^") + 27rilog(l - e^") 
we find that 

^(g) = ^(g) +7r. 

Since the shift by vr in the previous formula does not contribute to 9^i^, we conclude 
that ^J^L — ^HL. Finally, by working out several higher order terms can 
easily check that they always lead to the real functions 9, in accord with eqs.(5.24) 
and (5.28). 

Thus, we have shown that under the double Wick rotation the scalar factor 
remains unitary, while the dressing factor does not; the non-unitarity of the dressing 
factor is only due to the leading contribution 9afs, which is another distinguished 
property of 6'afs. Concluding this section, we note that it would be interesting to 
understand whether the BES factor [14] exhibits the same kind of non-unitarity 
behavior in the mirror theory. 
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6. Bethe ansatz equations 

In this section we discuss the nested Bethe equations for the hght-cone string theory 
on AdSs X and its mirror model. These equations based on the su(2|2) ©5u(2|2)- 
invariant string S-matrix [47] were recently derived by using the algebraic [60] and the 
coordinate [67] Bethe ansatz approaches. In the sector with even winding number, 
i.e. with the total momentum P satisfying e*^/^ = 1, the set of equations found 
in these papers coincides with the one previously obtained in [8, 9] by using the 
spin chain description of the gauge theory. It appears, however, that in the sector 
with odd winding number, where e^^^'^ = — 1, the Bethe equations by [60, 67] differ 
from the ones derived from the gauge theory. The origin of this disagreement can 
be traced back to the fact that in the light-cone gauge the fermions of the string 
sigma model are anti-periodic in the odd winding number sector [68, 56], and this 
changes the periodicity conditions for wave functions which one imposes to get the 
Bethe equations. Indeed, in the light-cone gauge one of the fields, an angle (p which 
parametrizes the five-sphere, appears to be unphysical and it is solved in terms of 
(transversal) physical fields. In particular, the equation of motion for (j) implies 



Since (p enters into parametrization of the five-sphere via e**^, the closed string peri- 
odicity condition for physical fields gives rise to the winding sectors 



where m is an integer. Now, we recall that fermions of the original string sigma- 
models are charged w.r.t. the U(l) isometry acting on as — + const. Also, 
the Wess-Zumino term in the sigma-model action contains e**^, i.e. it is non-local in 
terms of physical fields. To uncharge the fermions under the U(l) isometry, as well 
as to make the Wess-Zumino term local, one has to redefine the fermions as 



Thus, the redefined fermions acquire the periodicity properties which do depend on 
the winding sector 



i.e. they are periodic in the even winding sector and they are ant-periodic in the odd 
winding sector [68, 56]. 

As a result, the Bethe equations obtained in [60, 67] correctly describe the light- 
cone string theory in the sector with periodic fermions only. Changing the boundary 
conditions for fermions to anti-periodic, one derives a new set of Bethe equations 
which does agree with the gauge theory one for physical states satisfying e''^^'^ = — 1. 



(/)(27r) - 0(0) = P. 



4){2n) - 0(0) = 2nm, 
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6.1 BAE for a model with the su(2|2)-invariant S-matrix 

The asymptotic states of both the original and the mirror theory are constructed 
by applying the ZF operators to the vacuum state. The indices M and M 

are associated with two factors of the centrally-extended su(2|2) ©su(2|2) algebra; 
the latter being the symmetry algebra of the light-cone string theory [9, 10]. For our 
present purpose it is convenient to think about the ZF operator as being a product of 
two (anti) commuting operators each transforming in a fundamental representation of 
5u(2|2): A^^jjyj ~ ^M^]v;/- Since the string S-matrix is a tensor product of two su(2|2)- 
invariant S-matrices, the Bethe equations for the string model are, in a sense, the 
square of the Bethe equations for a model with the su(2|2)-invariant S-matrix. We 
start with discussing the Bethe equations for such a model. 

The multi-particle wave function which satisfies the Bethe periodicity conditions 
can be written as a superposition of the asymptotic states (see appendix 9.3.1 for 
details) 

\^) = J2'^^'---''-'AlM)---^liJPK'M^ (6.1) 

where is a number of particles in the asymptotic state and pi are their momenta. 
Denote by N{M) the number of particles of type M (that is number of indices of 
type M) occurring in the wave function (6.1). Obviously, 

= N{1) + N{2) + N{3) + N{A) . 

Since the scattering is elastic, the number of particles is a conserved quantity. 

The form of the Bethe equations derived through the nesting procedure of the 
coordinate Bethe ansatz depends on the choice of the initial reference state. Due to 
the su(2)^ bosonic symmetry there are two inequivalent choices for a model with the 
su(2|2)-invariant S-matrix. This is obviously related to the two forms of the Dynkin 
graph for su(2|2). 

First, one can choose a "bosonic" reference state which is created by acting with 
bosonic operators A\ on the vacuum: 

A\{p,) . . . Aiip^m ■ 

Then, we define 

K^^ = 2N{2) + iV(3) + A^(4) , i^'™ = A^(2) + iV(4) . 

It appears that in the scattering process not only but also these numbers are 
conserved [9]. By this reason, the values of and K^^^ are the same for any term 
in the sum (6.1). In particular, plays the role of the fermionic number, because 
in the background of the Aj-particles A\ counts for two fermions. The number K^^^ 
has a similar meaning. 
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Then the asymptotic Bethe equations based on the su(2|2)-invariant S-matrix 
for a sigma model on a circle of length R and with (anti)-periodic fermions can be 
written in the form [9, 60, 67] 





X7 T-r -Wi + 



Vk -wi-- 



(6.2) 



, , Wk — Vi + - -^-^ Wk — Wi — — 

1=1 g 1=1 9 

Ij^k 

Here e = for a sector with periodic fermions and e = 1 for a sector with anti-periodic 
fermions, xf depend on the momentum of the model, yi and wi are auxiliary roots 
of the second and third levels, respectively, and v = y + ^. 

On the other hand, if one chooses a "fermionic" reference state created by 
fermionic operators A^: 

Aipi) ■ ■ ■ AiPKm , 

then, one should define 

K^J: = 2iV(4) + N{1) + N{2) , i^'™ = N{2) + A^(4) , 



because these numbers are also conserved in the scattering process. Then, plays 

Ig-particles A\ 



the role of the bosonic number, because in the background of the Ag-particles A\ 



counts for two bosons. 

Then the corresponding Bethe equations take the following form 

^- + _ 

^v.R = {-lYl[So{pk,Pi)l[^^^\FT 

t\ t=i^k-yi^4 

l^k 

-ly = 1 1 n !:^z^ (6.3) 




l^k 

Equations (6.3) can be derived either by using the nesting procedure of the coordi- 
nate Bethe ansatz (see appendix 9.3.2 for an example) or by applying the duality 
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transformation discussed in [8] to eqs.(6.2). Comparing the two sets of Bethe equa- 
tions (6.2) and (6.3), we see that only the first hues in two sets are different. Let 
us stress, however, that in general ^ K^^. We further note that the bosonic 
reference state corresponds to, say, ?7i = 1 and the fermionic one to ?7i = — 1, where 
T]\ and ?72 are the gradings introduced in [8]. 

6.2 BAE based on the su(2|2) © su(2|2)-invariant string S-matrix 

The Bethe equations based on the su(2|2) © su(2|2)-invariant string S-matrix for 
both string and mirror models can be now easily written by taking a "product" 
of two copies of the Bethe equations for the su(2|2)-invariant model. Since any 
of the two sets, (6.2) and (6.3), can be used there are four different forms of the 
asymptotic Bethe equations based on the su(2|2) © su(2|2)-invariant S-matrix [8]. 
The corresponding bosonic reference states of the coordinate Bethe ansatz are of the 
form 



A;.(zi)...Aj.(z^i)|0), 7/1 = r/2 
and fermionic reference states are 



1; Aj.(^i)...4,(z^i)|0), 



33^ 



^71 = ^2 



-1 



13^ 



^1 



-ri2 



31^ 



^1 



-ri2 



where for the original theory the ^-variables lie on the real line, while for the mirror 
theory they have Im^; = ^2/21, and we also indicated the corresponding gradings. 

To discuss the bound states of the light-cone string sigma model, it is convenient 
to choose as the reference state the one created by the bosonic operators These 
reference states are dual to gauge theory operators from the su(2) sector. Then 
the corresponding Bethe equations based on the su(2|2) © 5u(2|2)-invariant string 
S-matrix can be written in the form [8, 9, 60, 67] 



SoiPk,Pi 



n 

1=1 

l^k 

(a) 

n Vk 
{a) 

1=1 Vk - ^. 



I ^k 
^k ^ 



^k 





rf. ' J. J. (a) (a) 



(a) 

n 



w 



(a) 



(a) 



K™ 



- V - 



1=1 ^k 



(a) 



n 



(a) (a) I 

^ ' vol' -w) + 



2i 
9 



g 1=1 ^k 

Ijtk 



(a) 



— W 



(a) 



2i 

g 



(6.4) 



Here we take into account that the string sigma model in the a = light-cone gauge 
is defined on a circle of length J, a = 1, 2 reflects the two copies of su(2|2) and yl°'^ 
and tW;"^ are auxiliary roots of the second and third levels, respectively, and v = y+-. 
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For the reader's convenience we point out that the excitation numbers in the set 
of Bethe equations are related to the ones used in [8] as follows 

and the Dynkin labels ^2] of su(4) and [si,S2] of su(2) © su(2) C su(2,2) are 
expressed in terms of the excitation numbers by the following formulas^^ 

q, = K'-Kll^, s, = Kll)-2Kl^, 

p=J-K'+ ^-{Kll^ + Kg)) , S2 = Kg) - 2Kg; , (6.5) 

To analyze the bound states of the mirror theory, it is more convenient, however, 
to choose as an initial reference state the one created by the operators A^^^ . The 
reason is that the operators create states from the sl{2) sector, and, as we have 
seen, it is this sector which gives rise to mirror magnons. Analogously, there are 
M-particle bound states made only out of the A^^-type particles. 

If we choose in the mirror theory the above-described reference state then the 
corresponding Bethe equations take the form 

2 ^fa) 

e^^~'=^=n[^o(p.,pof nn 

1=1 a=l 1=1 





(6.6) 



(a) (a) 2i ' 

Note that in the mirror model we do not have (—1)*^ in the middle equation because 
the fermions are always anti-periodic^^ with respect to a. In terms of excitation 

^^Lct us note in passing that in recent papers [70, 71] the anomalous dimension of the operator 
TrJ-^ was computed by using the asymptotic Bethe ansatz with an understanding that in the large 
L limit one may trust the corresponding result to an arbitrary loop order. One can notice, however, 
that the excitation pattern of Bethe roots for the operator is {Kjl^ , i^Tg^j , , iCgj , A"gj) = 
(0,2L - 3,2L - 2,2i - 4,i - 2) with J = f , and, therefore, one would expect the breakdown 
of the asymptotic ansatz due to the finite size effects already at two loops. It may happen that 
the asymptotic ansatz could still be used to determine the leading L behavior of the anomalous 
dimension of TrT^ if the finite-size corrections are subleading at large L, but this is currently 
unknown. 

^^We are grateful to R. Janik and M. Martins for drawing our attention to this point. 
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numbers, the Dynkin labels read now as follows 

p=J- ^-{Kll^ + Kll)) + K}^ + K}^ , s, = K'- Kll) , (6.7) 

12 - J^(2) - -"-(2) • 

7. Bound states of the AdSs x gauge-fixed model 

Bound states arise as poles of the multi-particle S-matrix corresponding to complex 
values of the particle momenta, see e.g. [69]. In the thermodynamic limit they 
are described by string-like solutions known as "Bethe strings". In this section we 
discuss in detail the bound states of the string sigma-model. They have been already 
analyzed in [50, 51]. The main outcome of this analysis is that the M-particle bound 
states comprise into short (BPS) multiplets of the centrally extended su(2|2)©su(2|2) 
symmetry algebra. Although the S-matrix exhibits additional simple and double 
poles beyond those corresponding to the BPS multiplets, these singularities do not 
lead however to the appearance of new bound states [52]. In other words, the only 
bound states in the theory are the BPS ones. As we will see, they exist for all 
values of the (real) bound state momentum — tt < p < it, but have a rather intricate 
structure. Moreover, depending on the choice of the physical region for a given value 
of the bound state momentum there could be 1, 2 or 2^^"^ M-particle bound states 
sharing the same set of global conserved charges: Qr = YliLi Qri^i) ■ It is unclear 
to us whether this indicates that the actual physical region is the one that contains 
only a single M-particle bound state (it is the one with \x^\ > 1) or it is a sign of a 
hidden symmetry of the model responsible for the degeneracy of the spectrum. 

7.1 Two-particle bound states 

Let us consider a bound state made of two excitations from the su(2) sector of the 
string sigma-model. In terms of the ZF creation operators we can think about this 
state as 

where the particle momenta pi and p2 are complex. We find it convenient to 
parametrize the momenta as follows 

p p 
Pi = ^-i(l, P2 = ^ + iq, Reg>0, (7.1) 

where p is the real total momentum of the bound state. When q is real then pi and 
P2 are complex conjugate to each other and the energy of the corresponding bound 
state being the sum of the (complex) energies of individual particles is obviously real. 
Interestingly, as we will show below, there necessarily exists a branch of BPS bound 
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states which corresponds to complex values of q with Re g > 0. Such solutions can 
be reinterpreted as solutions parametrized by a new real variable q: q —>■ Keq and 
for which the real parts of pi and p2 are not anymore equal to each other. Of course, 
one has to check that the energy of these solutions is real. 

The first equation in the set of the Bethe equations [6] takes the form 

+_ _1_ 1 

g^(p/2+Im,)LgRegL ^ ^iP TT '-1^ , (7.2) 

1=2 1 ' x+ 

where P = pi + p2 + ■ ■ ■ +Pk^ and L = J + with J being one of the global charges 
corresponding to the isometries of the five-sphere. 

We see that for large L the l.h.s. is exponentially divergent. Then, there should 
exist a root p2 such that for Re g > we have^^ 

In the infinite L limit eq.(7.3) becomes 

(x^-xt) (l-^) =0, (7.4) 
V X1X2 ^ 

which is equivalent to 



-xt = or 1 = . (7.5) 



Xi X2 

The first equation 

x--xt = (7.6) 

implies that the central charges corresponding to the two-particle bound state satu- 
rate the BPS condition [50] 

i72 = 22 + 4/sin2|. (7.7) 

On the contrary, solutions of the second equation in (7.5) do not saturate the BPS 
bound, and as was argued in [52], this pole of the S-matrix does not correspond to a 
bound state. 

It is easy to see that equation (7.6) is equivalent to vanishing the following fourth 
order polynomial^® in the variable t = cos | 

Ae\t - e'?)(l - eH) + - 1)(1 - 2eH + e^'^f = . (7.8) 

^^We assume here and in what follows that the dressing factor 1712 is non-singular on solutions of 
the bound state equation. 

^^For any p there are two solutions for and, therefore, for x+ = e^^x^. The fourth order 
polynomial is universal and it does not depend on which solution for x^ we take. 
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The equation has four solutions which can be cast to the following simple form 

[sjl + g"" sin^l ± [cos I sjl + g"" sin^ £ ± ^cos^ £ - g^ sin^ f) 

g"^ siv? £ ^) 

where any choice of the ± sign is possible. 

Analysis of eq.(7.9) immediately shows that solutions corresponding to real val- 
ues of q exist if and only if the total momentum p does not exceed a critical value 
Per determined by 



sin^ = - l) • (7.10) 



2 Per _ 



For any given p obeying \p\ < pcr equation (7.8) has four real roots for q, two of them 
are positive and the other two are negative. According to our assumption Re g > 0, 
only positive roots are acceptable. They are given by the formula 

{^1+9^ sin^l + l) (cos £ yr+v"^ ± y/cos2 £ - g^ sin^ 

sin^ £ ^ 

Various expansions of eq.(7.11) for small and large values of g can be found in Ap- 
pendix 9.4. 

It turns out that from the two positive roots only the smaller one, g_, falls inside 
the region confined by the curves \x^\ = 1. We therefore arrive at the conclusion 
that inside the region > 1 there is a unique solution with real p and q, and it 
exists if and only if 



20 



2g+ \ 2W1+Ag^ ~2 

bl < Per , < g < log y . (7.12) 

v/l + 4^2 _ 1 

The second solution with q = q^ lies outside the region with |a;^| > 1 but inside 
the region with —002/21 < Im(2;) < 002/21; the latter maps onto the complex p-plane, 
see section 4. Both solutions have the same values of all global conserved charges 
Qr = qr{zi) + qr{z2) = ^ [(^i")^"*" ~ (2^2)^"''] bccausc and X2 are the same on 
both solutions. 

We see that if we choose the physical region to be the one with > 1 then 
there is a unique bound state with \p\ < p^r- This region, however, does not cover 
the whole complex p-plane. One the other hand, if the physical region is the half of 
the torus corresponding to the p-plane, then there are two solutions with the same 

^'^Thc solutions with negative q correspond to bound states of anti-particles with negative energy. 
^"The energy of the corresponding bound state is < \J2^/Y+ Ag^ + 2 . 
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energy and other conserved charges. Finally, if one considers solutions on the z-torus 
then there are four solutions but only two of them have positive energy. 

Continuing above the critical value, \p\ > pcr, two solutions (7.11) acquire imag- 
inary parts and become complex-conjugate to each other, or, equivalently, the real 
parts of pi and p2 become different. Thus, we see that the BPS bound states naturally 
split into two families depending on whether the total momentum is below (the first 
family) or above (the second family) the critical value pcr- 

The two complex conjugate roots give two different solutions beyond criticality: 
pf = -±lmq — iReq , pf = -^Imq + iReq , Reg>0. (7.13) 

We can choose either {pt,pt) or (p]~,p^) as a possible solution of the BPS condition 
(7.6). Note that the second solution is the complex conjugate of the first one. A re- 
markable fact to be proven below is that both solutions lie precisely on the boundary 
of the region defined by the curves | | = 1 . 

Now if we adopt the physical region (sheet) to be |a;^| > 1 with the boundary 
|a;^| = 1, then it should contain only one solution from the second BPS family. 
Indeed, we do not expect the doubling of the number of BPS bound states moving 
beyond the critical point. The second solution can be then naturally interpreted as 
lying on the boundary of another (unphysical) sheet joint to the physical one through 
the cut. It is unclear however what is the precise origin for such an asymmetry. A 
possible explanation would be the absence of parity invariance of the string sigma- 
model but a concrete implication of this fact is unknown to us. 

To visualize the singularities of the string S-matrix and also to verify that energy 
is real for the second BPS family, it is instructive to analyze eqs.(7.5) in terms of 
the generalized rapidity variables zi and Z2 associated to the first and the second 
particles, respectively. It is not hard to see that the first family of the BPS states 
corresponds to imposing the reality condition Z2 = zi. In this case, eqs.(7.5) are 
equivalent to 

Im(x^) = or |x|f| = l, (7.14) 

where the first equation defines the first BPS family. Solving the bound state equa- 
tion for Zi, one gets a curve in the torus. The part of the curve that lies inside 
the region > 1 is represented in Fig. 5 by the green curve BiOCi, and the 
corresponding momentum pi has Im(pi) = — g_. The variable Z2 = z\ of the sec- 
ond particle runs along another (conjugate) green curve B2OC2, which can be also 
viewed as describing solutions of the equation Im(x^) = for Z2. The dashed curves 
on Fig.Sa, which are outside the region > 1, represent solutions of the equa- 
tions \m{x^) = Im(a;^) = for Zi,Z2 corresponding to the momentum pi with 
Im(pi) = 
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a) b) c) 



Figure 5: Two-particle bound states of string theory. Figure a) describes the first BPS 
family corresponding to p < pcr- The green curves are Im(x~) = for lm{z) < and 
Im(x"'") = for Im{z) > 0. For any p < p^r there are two solutions: the first one is 
represented by the continuous curves BiOCi (1st particle) and B2OC2 (2nd particle), the 
second one corresponds to the dashed curves AiBiUCiDi (1st particle) and A2B2UC2D2 
(2nd particle). Figure b) describes the second BPS family corresponding to p > pcr- 
Again, for any p > pcr there are two solutions B2C2 U AiBi U CiDi and BiCi U A2B2 U 
C2D2. Figure c) represents one of the four possibilities to smoothly connect solutions from 
the first and the second BPS families. Here the variable zi of the 1st particle is on the 
curve AiBiOCiDi. When zi runs along the curve from Ai to Di the real part of the 
momentum Re(pi) increases monotonically from — vr to it. At the same time, the variable 
Z2 corresponding to the 2nd particle encloses the curve A2B2OC2D2. 

To describe the second family of the BPS states corresponding to the complex 
values of q one has to take 

^2 = -zl + ^ + ^. (7.15) 



2 2 



In this case 



where we have used the properties of Jacobi elliptic functions under the shifts by 
quarter-periods. Hence, due to the BPS equation = x^, the points zi and Z2 lie 
on the curves = 1 and \x'^\ = 1, respectively. 

As was discussed above, there are two different ways to choose the second BPS 
family which is equivalent to deciding what is the physical sheet. Consider the point 
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Zi corresponding to the first particle, Fig.Sc. When it moves along the curve BiOCi 
corresponding to the first BPS family and reaches, e.g., the point Ci then there are 
two possibilities to continue its path along the curve |a;~| = 1: either one moves 
along CiDi or along CiBi. In the case when zi moves along the curve CiDi, the 
second point Z2 follows the path C2D2. In the opposite situation, when zi moves 
along CiBi, the point Z2 follows C2B2. Similar discussion applies to continuing the 
first BPS family beyond Bi. Obviously, for the second family zi and 22 are not 
complex conjugate anymore, rather they obey the relation (7.15). The bound state 
energy H = ig{x2 — x^*") — 2 is however real, as one can also check by using the 
shift/reflection properties of the elliptic functions. 

Our discussion reveals that there are four special points on the 2;-plane 

Zcr = ±—±— (7.17) 

where both equations Im(x]") = and \xi \ = 1 or Im(x^) = and = 1 are 
simultaneously satisfied. These are the critical points where two BPS families meet. 

The most transparent description of the bound states is achieved in terms of the 
rapidity variable u introduced in section 4, rather than in terms of momentum p or 
the variable z. Indeed, in terms of u eq.(7.3) becomes 



{x^ -xt) (1- ^-^) =ui-U2-- = 0, (7.18) 
\ Xi xj / g 



i.e. the rapidity variables Ui and U2 of the first and the second particle lie on a 
straight line running parallel to the imaginary axis. Moreover, for the first BPS 
family the variables ui^2 are subject to the following conjugation rule ul = U2 which, 
together with eq.(7.18) allows one to conclude that 

Ml2=Mo±-, MoeM. (7.19) 

9 

This is a typical pattern of "Bethe string". One can further see that for the first 
BPS family corresponding to p < p^r the variable Uq is restricted to satisfy 

|mo|>2, Micr = ±2+-, (7.20) 

9 

where ui^cr is a critical value of rapidity ui for which the first BPS family ceased to 
exist. Under the map to the w-plane the four critical points mapped to the 

four branch points on the M-plane (see Fig. 3 in section 4) 

Mcr = ±2 ± - . (7.21) 
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Let us now turn to the second BPS family. First, by using eq.(4.5) and the 
properties of the elhptic functions, we derive 



«i = X- ( - ^2 + ^ + ^) + —f ^ ^ + - = W2 . (7.22) 

^ 2 2 / (■ _ ^2 + f + f j ^ 

We see that for both famihes of BPS states the conjugation rule for m's is the one 
and the same. By this reason, a solution to the BPS condition is always represented 
by the Bethe string (7.19). However, one finds that for the second family a solution 
exists for \uq\ < 2 only. Thus, on the w-plane both families of BPS states admit a 
uniform description in terms of the Bethe string with mq running over the whole real 
line. 



7.2 Multi-particle bound states 

The consideration of the two-particle bound states can be easily generalized to the 
M-particle case. The corresponding set of bound state equations reads [50] 

x--x++i=0, j = l,...,M. (7.23) 

The total momentum of a state satisfying these equations is given by 

,T^H~ /->-»~l~ /Trl I 

ip _ -^-^ -^M _ -^1 

6 _ _ ■ ■ ■ _ _ ; 

and the energy of the state is 

M 

if^/ = ^ (-1 - igx^ + igxl) = -M - igxf + igx^j . (7.24) 

Both the energy and momentum depend on the values of xf and x]^,j only. Since 
the energy is real, a;^^ must be the complex conjugate of xf: {x~^)* = xf. In fact, 
a simple but important observation is that any global conserved charge of a state 
obeying (7.23) depends only on xf and x^: 

MM. 

i=l i=l 



Another important consequence of eqs.(7.23) is that the coordinates x^ and x^j 
satisfy the following quadratic constraint 

+ 1 - 1 2M. 

xj + — - Xj^j — = 1 . (7.25) 

Xi 9 
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This is the same constraint as the one satisfied by x^ (4.3) with g —>■ g/M, and we 
get immediately the dependence of xf and xjj on the total real momentum p 



6 2 



X 



X 



2g sin - 



M 



2g sin 2 

and, using (7.24), the BPS energy formula 



M2 + 4^(2 sin' 



P 



(7.26) 



H 



,, = M2 + 4/sin2| 



Moreover, we see that the set of global conserved charges Qr is the same for any 
solution of (7.23) with a given total momentum p. 

It is also easy to see that the number of different solutions with a real momentum 
p and positive energy is equal to 2^^~^ because for a given x^ there are two different 
x~ solving the constraint (4.3), see the diagram below for M = 4 



X, 



X 1 



X n 



X n 



X n 



X n 



X. 



Xn 



X q 



Xn 



X3 


— 3^4 




— 3^4 


^3 


= 


^3 


= X^ 


^3 


= X^ 


X3 


= X^ 


X3 


= X4 


X3 


= X4 



Xa 



To have all these solutions one would have to allow the parameters Zi of the particles 
to be anywhere on the z-torus, in particular, some of them would be in the anti- 
particle region with |x^| < 1. 

However, if we require that all the constituent particles of the bound state belong 
to the region |x^| > 1 then we are left with a unique solution because for a given 
x^ only one solution for x~ satisfies the condition |x~| > 1. For M even it is also 
necessary to specify what parts of the boundaries |x^| = 1 belong to the region 
because if the momentum of a bound state exceeds a critical, g— and M-dependent, 
value then there are several solutions of the bound state equations with jx^^^gl — 



X 



M/2+1 1 



1. 



Finally, if the parameters Zi of the particles belong to the half of the torus 
corresponding to the complex p-plane, then one can show that for any M there are 
two solutions of the bound state equations. 



^^In general for a given momentum p there are two solutions of the constraint (7.25), and there 
could be any sign in front of the square root in (7.26). The positive sign guarantees the positivity 
of the energy. 



-52- 



Just as for the case of two-particle bound states, the simplest description of M- 
particle bound states is provided by the w-plane where a solution is given by the 
Bethe string 



= uo + {M- 2j + 1) 



j = l,...,M. 



(7.27) 



We can choose one and the same map of the u-plane with the cuts described in 
section 4 onto the region of the z-torus with > 1 for all the particles. It is 
then obvious that for a given momentum p there is just a single M-particle bound 
state that falls inside the physical region. Its structural description however becomes 
rather involved. 

7.3 Finite-size corrections to the bound states 

It is of interest to analyze finite-size corrections to the energy of the BPS bound 
states, and to see what restrictions on the dressing factor could be derived from the 
condition that the energy corrections are real. To this end, we consider two-particle 
states in the 5u(2) sector described by the following two equations, see (7.2) 



where S12 = ^+^- ai2 is the unitary factor that appeared in the splitting (5.25) of the 



scalar factor, and J is one of the global charges corresponding to the isometries of 
the five-sphere. The variables xf also satisfy the constraint (4.3). These equations 
are supposed to be valid asymptotically for large values of J, and have to be modified 
for finite J. 

We will analyze these equations for large values of J in the vicinity of a bound 
state satisfying the bound state equation = and having a fixed total mo- 
mentum p = where m is an integer. The quantization condition for the total 
momentum follows from the second equation in (7.28). 

Let xf denote the values of xf satisfying the bound state equation and the 





.+ 



second equation in (7.28). Then, 




exponentially decreases at large J, 



and we can look for a solution of the form 




Expanding the equations (7.28) and the constraint (4.3) in powers of , we find a 



system of linear equations for yf. The solution of the system is given below 
AE,2{tg{^2 - x^) + X2 x+) (2^x+X2+ + g (x+ ((x+)2 + 1) - 2X2+)) 



yt 
yt 



(72(^X^-(^ + 2^X2-)x+) ((x+)2-l)' 
4Si2X^('i5((x2 - x^) + X2 x^) 



^7((7X2--(^7 + 2zx2)x+)((x+)2-l) 
4iSi2X2 (5'(X]^ - x^) + iX]^x^) 



+ 22X2 )X+- (7X2 )((x+)2-l) 

4Si2(^(xj^ - x^) + ixj^x^) (^X2 (x^)^ - 2{g + 2X2 )x^ + g:>q) 
g^-igi:^, - x+) - 2x2 x+) ((x+)2 - 1)' 



where E12 is evaluated on the solution to the bound state equation. The leading 
correction to the energy of the state is easily found by expanding 

E = Ei + E2, Ei = l + ^ - ^ = -l-igxf + igx; . (7.29) 
By using Ei = 1 + ^ - we obtain 

/x7\^ 42(x^(2x+ -x+) -X+X+) 

VX+ y (^(X2 - X+) - 2^X2 X+) ((x+)2 - 1) 

On the other hand by using Ei = —1 — igxf + igx~ ^ we get 

4^X2 X]^X^(X2 + x+ — 2x^) 

^ = U J ''^^,(x2--X+)-2zX2-X+)((x+)^-l) • ^'-^'^ 

Even though the expressions look different they coincide on solutions to the bound 
state equation. In what follows we will be using the simpler eq.(7.31). Note also that 
the perturbation theory breaks down at p = pcr- Due to the quantization condition 
for the momentum p it may happen only at special values of the coupling constant 
g depending on m/J. 

It is clear that the energy correction cannot be real for any choice of the dressing 
factor S12. The imaginary part of the correction depends on the branch of the bound 
state under consideration. 

In the first case with Im(x^) = Im(x^) = and the total momentum smaller 
than the critical value (7.10), the parameters x+ satisfy the complex conjugation rule 




5E-6E''=\{'-^\ T.,,-{'^\\\} 42X2-x+x+(x2- + x+ 2x+ 



(^7(X2 -X+)-2zX2X+)((x+)2-l) • 
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Taking into account that (^~f^ ~(^) conclude that in this case the correc- 
tion is real only if the dressing factor is real ZI12 = S*2 . This property of the dressing 
factor can be easily shown by using the representation (5.22) for the dressing phase. 

In the second case with \x.]^\ = |x^| = 1 and the total momentum exceeding 
the critical value (7.10), the parameters satisfy the complex conjugation rule 
(x^)* = X2 , (xj")* = 1/x^, and we obtain 

6E - 6E* = 

_ 4X2-X+ (El2(x+X2-)-^(2 - (x^- + x+)x+) - Si2(x+)-^(x+)-^+Hx2- + x+ - 2x+)) 

(Z^7(X2 -X+) + 2X2X+)((x+)2-l) 

We see that the imaginary part of the correction would vanish only if 

2xi 



^12 ~ Sl2(x^ 



\2J+1 V^2 "T ^1 ^^2 



(2 



Since the last equation depends on J and on a particular bound state solution, it 
cannot be satisfied for any choice of the dressing factor. The complex energy of the 
state would mean that the Hamiltonian of the model is not hermitian for finite J. 

One might conclude from this result that the S-matrix poles with |x^| = |x^| = 1 
are spurious and do not correspond to bound states, and, therefore, should be omit- 
ted. That would mean, however, that for any value of the total momentum the bound 
states satisfying the equations Im(x^) = Im(x^) = would disappear as soon as the 
coupling constant g reaches a critical (momentum-dependent) value. This seems to 
contradict to the statement that the bound states are BPS. We believe that such 
a conclusion might be erroneous and the result above indicates, in fact, that the 
asymptotic Bethe ansatz cannot be used to analyze the finite-size corrections to the 
energy of bound states with the total momentum exceeding the critical value (7.10). 

To show that this is indeed the case, let us recall that, as was shown in [38], 
at large values of the string tension g and the charge J the dispersion relation re- 
ceives finite-size corrections of the order e~'^^^^^^^^^'^\ On the other hand, the energy 
correction we computed above is of the order e""^"^ where q is the imaginary part 
of the momentum p2. It depends on the total momentum p and the string tension 
g. By using eq.(9.35), it is not difficult to determine the large g dependence of the 
momenta pi and p2 of a bound state 

cos f i I s cos f i / 1 X / X 

= o 2 ■ 3p - — ^ + ^(-^) ' P2=P- ^ , .\r> + — ^ + ■ (7-32) 

2g^ sm I sm I g-^ 2g^ sm | sm | g-^ 

The second solution of eq.(7.8) (with g > 0) is related to (7.32) as Pi P2 , P2 ^ 
pI that is one exchanges the real parts of momenta pi. A surprising result of the 
computation is that q is equal to ^ p , and, therefore, e"'^"' is exactly equal to the 
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magnitude of the finite-size correction to the dispersion relation. That means that 
computing the finite J correction to the energy of such a bound state one has to 
take into account the necessary modifications of the asymptotic Bethe ansatz. As a 
result of these modifications, one should be able to get a real finite-size correction 
to the energy of a bound state carrying momentum exceeding the critical value. In 
fact, this would be a non-trivial check of finite J "Bethe" equations. 

The analysis performed above raises the question if one can use the asymptotic 
Bethe ansatz to compute the corrections to the energy of the bound states with 
momenta smaller than the critical value. At large g we can again compare the value 
of q with p . If q is less than ^.^^ p then the energy correction (7.31) is bigger than 
the correction due to finite J modifications of the asymptotic Bethe ansatz, and we 
can trust (7.31). Since pcr = at large values of g one should consider a bound 

state with momentum p of the order l/^/^- The large g dependence of the momenta 
Pi and p2 of a bound state is easily found by using eq.(9.36) 



p. = l-2 '^^'- " +.... + " +.... (7,33) 

2 gp 2 gp 

leading for p < pcr to the following two real solutions for q 



1 ± \/l — 2^ 

= 2 ^ — . (7.34) 

9V 

Comparing these values with — ^ ^ — , we see that (/_ < — and (/+>—. Thus, the 
asymptotic Bethe ansatz can be used to analyze finite J corrections to the energy of 
a bound state with momentum smaller than p^r for the bound state with q_ only. 

Actually, the fact that the energy correction (7.31) to the bound state with qj^ 
is smaller than the corrections due to finite J modifications of the asymptotic Bethe 
ansatz raises a question if these solutions correspond to the actual bound states. It 
may happen that finite J Bethe equations would not have any solution that would 
reduce to the solution with g+ in the limit J — > 00. 

A similar analysis can also be performed for small values of g. Then we expect 
that the finite J effects (in gauge theory they are due to the wrapping interactions) 
become important at order g"^^ , and therefore we could trust the asymptotic Bethe 
ansatz and the energy correction (7.31) only if g < — 21og5'. 

The leading small g dependence of q of the bound state solutions with the mo- 
mentum smaller than p^r is given by eqs.(9.37), (9.38) 

p 

q+ = ~2\ogg^ , g_ = -logcos-H . (7.35) 

We see immediately that again only the solution with the smaller imaginary part of 
the momentum g_ satisfies the necessary condition. The energy correction to the 
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state with g+ is of order g"^ that is exactly the order of wrapping interactions, and 
the asymptotic Bethe ansatz again breaks down for the state. 

Finally, the leading small g dependence of q of the bound state solutions with 
the momentum exceeding p^r is given by (9.39) 

g± = -log| ± m H , (7.36) 

where a is related to the momentum p as follows p = tt — 2g cos a . 

We see that the real part of q± is smaller than — 21og5f, and therefore one could 
conclude that one might use the asymptotic Bethe ansatz for the states in this regime. 
This, however, leads to the problem of the complex energy of these states discussed 
above. As before the only resolution of the problem we see is the breakdown of 
the asymptotic Bethe ansatz. This would imply, however, that for these states the 
wrapping interactions become important already at the order g'^ . The fact that in 
gauge theory these states are not dual to gauge-invariant operators does not seem 
to be important for this conclusion. One could for example scatter a bound state 
carrying momentum p = n which always exceeds the critical momentum pcr with an 
elementary one carrying momentum — tt so that the total momentum would be zero, 
and such a state would be dual to a gauge-invariant operator. We would still expect 
the finite J corrections to this state to be of the order g"^ . Another puzzling property 
of the states with p > p^r is that in the limit g ^ the states are pushed away from 
the spectrum because pcr = ti", and cannot be seen in the perturbative gauge theory. 



8. Bound states of the mirror model 

Let us now consider in a similar fashion bound states of the mirror model. In this 
case one should consider mirror particles of type A^^. 

We begin our consideration with two-particle bound states, and let the complex 
momenta of the two particles be 

Pi = ^-iq. P2 = ^ + iq, Re g > , 

where p is the total momentum of the mirror bound state. 
The first equation in (6.6) takes the form 

„+ 1 - ^ 

^2 

Xn 1 — 



where we set all auxiliary roots to 0. Assuming that the dressing factor does not 
vanish, we conclude that for Re g > and in the limit R ^ oo the following bound 
state equation should hold 

xt-Xo=0. (8.2) 
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The second factor in the denominator of the Bethe equation (8.1) may also vanish 
but the energy of the corresponding state does not satisfy the BPS condition. We 
expect that, just as a similar factor in the string theory, the pole due to this factor 
does not correspond to a bound state. 

By using eqs. (2.15) which express as functions of we find that eq.(8.2) is 
equivalent to 

- + - 2qH (2 - t) + g^t = , (8.3) 

2 

where t = 1 + This equation gives the following two solutions with a positive real 
part of q: 

Solutions for q are real provided the expression under the square root is nonnegative, 
and this implies the following restriction on the total momentum of the bound state 

\p\ < Per = V2^-1 + Vl+ V • (8.5) 

For an exact inequality we have two positive solutions q±, and when the bound on 
the momentum is saturated the solution is obviously unique^^ 

g_ < gcr < q+ , Qcr = -^^1 + V_ (g g) 

It is interesting to notice that the dependence of q± on the momentum of the bound 
state is smoother at p = than the one for string theory bound states. We see from 
eq.(8.4) that g_ reaches its minimum, and g+ reaches its maximum at p = 

g-- = yrT7-^7, gr' = A/r+7 + ^7, p = o. (8.?) 

In string theory the corresponding values are and oo. 

To find what curves in the z-torus correspond to the two solutions with real q 
we take into account that in this case pi* = P2, and the reality condition for in 
the mirror theory is {x^Y = l/x2 - Thus the bound state equation (8.2) reduces to 
the following equivalent conditions 

\x~l\ = 1 < > \X2 1 = 1, 

being represented by the two curves in the 2-torus that bound the yellow region 
with Ix'^l < 1, > 1 in Figure 1. Note that the curves are symmetric about 

^^The energy of the bound state is £„■ = 2 arcsinh^ \Jl + ly^l + ig^ . 
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the horizontal hne passing through the point z = Let us recall that hermitian 
conjugation in the mirror theory is defined with respect to this line, see section 4. 
It is not difficult to check that the parts of the curves \xi\ = 1 , | = 1 that are 
inside the region Im(x^) < correspond to the smaller root g_ of eq. (8.2). The 
other parts of the curves correspond to the second solution with g = g+, see Figure 
1. Just as it was for string theory bound states, both solutions have the same values 
of all global conserved charges Qr = qr{zi) + qr{z2) = ^ [ — {x^Y"^ + (a;^)^"''] . 

We see that if we want to have only one bound state with \p\ < p^r in a physical 
region, then we should choose the physical region to be the one with Im(a;^) < but 
not the one bounded by the curves = 1 as it is for string theory. We will see 
in a moment that the region Im(x^) < also contains bound states with \p\ > pcr 
described by the solutions with complex q. 

Above the critical value, \p\ > pcr, the two solutions (8.4) acquire imaginary 
parts and become complex conjugate to each other. It is convenient to denote the 
corresponding solutions as follows 



q± 




p2(4 + p2) 



We see that the real part of q± is a decreasing function of p, and its minimum value 
is 1. On the contrary the imaginary part of q± is an increasing function of p and it 
behaves as ±p/2 at large values of p. As a result, the two complex momenta 

Pi^ = -±Img — iReg, p2^ = - ^ Im q + iRe q , Reg>0, 

have the following large p behavior 

P\^ = p - i, p2^ = i ; pi^ = -i , p2~ = p + i. 

A remarkable fact is that both solutions lie precisely on the boundary of the 
region Im(x^) < 0. To see this we notice that, just as it was for string theory bound 
states, the coordinates Zi and Z2 of the solutions with the complex values of q are 
related by eq.(7.15) 

^2 = -^1* + y + y . (8.9) 
Then, one can easily show that 

x-{z2) = ^" ( - ^1* + y + y) = ^"(^n = [^^(^i)]* ' 

and, therefore, the bound state equation xf = X2 is equivalent to Im(a;"''(2;i)) = 
Im(x~(2;2)) = 0. We plot the corresponding curves in Figure 6. 
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a) b) c) 



Figure 6: Bound states of the mirror theory. Figure a) represents the first BPS family: for 
any p with \p\ < pcr there are two solutions corresponding to (the curves BiCi for the 1st 
particle and B2C2 for the 2nd one, respectively) and to g+ (the curves AiBiUCiDi for the 
1st particle and A2B2UC2D2 for the 2nd one, respectively). Figure b) represents the second 
BPS family which is also doubly degenerate: it is given by either AiBi U CiDi U B2D2 or 
by BiCi U A1B2 U D1C2. Figure c) corresponds to one of the four possibilities to connect 
the first and the second BPS family: when the variable zi of the 1st particle runs along 
the curve AiBiCiDi the real part of its momentum increases from —00 to +00. At the 
same time, the variable Z2 of the 2nd particle encloses the curve A2B2C2D2. 



Thus, we have shown that these solutions lie on the boundary of the region 
Im(x^) < 0, and, therefore, the region contains bound states with any value of the 
total momentum and could be considered as the physical one for the mirror model. 
It is also necessary to specify what part of the boundary of the region Im(a;^) < 
belongs to the physical region, and this can be done by choosing properly the cuts 
in the u-plane where the bound state equation reduces to 

, _ 2i 
xT — Xo =0 =^ U2 — ui = — . 

9 

As was discussed in section 4, eqs. \xf\ = | = 1 describing a bound state with 
the momentum not exceeding the critical value Pcr and with a real q give a Bethe 
string solution with the real part of u lying in the interval [—2,2] 

i 

ui,2 = uoT - , -2<Mo<2. 
9 

On the other hand, values of uq lying outside the interval [—2, 2] correspond to 
solutions of eqs. lm.{xi) = Im(x2 ) = 0. The momentum p = p{u) is a multi-valued 
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function of u, and one should choose a proper branch of the function to get the right 
values of the momenta pi , p2 of the bound state. This fixes the cuts in the w-plane 
which run from ±00 to ±2 =p -, and also the boundaries of the region Im(x^) < m 
the z-plane which is mapped onto the w-plane with these cuts. 

The discussion of bound states of M particles of type A^^^ basically repeats the 
one in section 7. One finds a system of equations 

x+-xT^, = 0, j = l,...,M-l. (8.10) 

In terms of the variable u the Bethe string solution reads as 

u, = no - (M - 2j + 1)^ , J = 1, . . . , M , (8.11) 

and has the energy 

£" = log ^ = 2 arcsinh ^ + p2 ^ (8.12) 

where p = pi + . . . + pm is a total (real) momentum of the bound state. 

Depending on a choice of the physical region, the system (8.10) could have one, 
two or 2^^"^ solutions. All solutions have the same global conserved charged. They 
behave, however, differently for very large but finite values of R, and the solutions 
which are not in the region Im(x^) < show various signs of pathological behavior. 
In particular, they might have complex finite R correction to the energy, or the 
correction would exceed the correction due to finite R modifications of the Bethe 
equations thus making the asymptotic Bethe ansatz inapplicable. 
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9. Appendices 



9.1 Gauge-fixed Lagrangian. 

The Lagrangian density of the gauge-fixed sigma-model in the generahzed a-gauge 
[56, 57] can be written in the following form [26] 



a Gtt + (1 — 0.) G. 



ipip 



1 - aYG^^ - a^Gtt (1 - a^G^^ - a^Gu 



(9.1) 



where 



W = 1 



[1 — aY'Ginip — ci^Gtt 



G^pipGtt 



+ 



GiptpGtt 

((1 - afG^^ - a'Gttf 
2GiptpGtt 



X ■ X + X' ■ X' ^ 

{{d^X ■ d'^Xy - {d^X ■ dpXf) 



Here X = {y'^,z^), where y*, i = 1,...,4 are four fields parametrizing five-sphere, 
while are fields parametrizing four directions in AdSs. The fields X in the La- 
grangian above are contracted with the help of the metric 

ds'^ = -Gttdt^ + Gzzdz^ + G^^dif'^ + Gyydy^ . 



Here 



Gtt 



l + z' 
l-z^ 



2\ 2 



(1-.2) 



2 ' 



Gu 



l-y 



2\ 2 



1 + r 



2 ] ' ^yy 



(l + y2) 



2 ' 



where we had used the notation z = z'^z^ and y = y^y\ 



9.2 One-loop S-matrix 

Here we describe the properties of the "one-loop" S-matrix which is obtained from 
the S-matrix (3.24) upon taking the limit g ^ 0. We will work in the elliptic 
parametrization discussed in section 4.1. According to eq.(4.18), in this limit Jacobi 
elliptic functions degenerate into the corresponding trigonometric ones and we find 
the following trigonometric S-matrix: 
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S{zi,Z2) 



-i{zi-Z2) 



cot zi — cot Z2 + 2i 



cot Z\ — 


cot Z2 — 


Zii 


i(zi-Z2) 


2i 




cot Zi — 


cot Z2 — 


On 






— 


cot Zi — 


cot Z2 — 


2i 


cot 2:1 


— cot Z2 




cot Zi — 


cot Z2 — 


2i 




— cot Z2 




cot zi — 


cot Z2 — 


2i 


1(21-22) 


2i 




cot Zi — 


cot Z2 — 


2i 


1(21-22) 


2i 




cot zi — 


cot Z2 — 


2i 



El (g) El + E 
El (E)E^ + E. 



EI(^Ei 



E. 



El + El ®El + El® El 

El - El (g) El - El (g) El 
E^ 



El 



' E^ + E^ 



E: 



E^®El + El®El~ El 
El ®El + El ®El + El®El + El 
El (g El + El® El + El ®El + Et®E. 
El g) El + eI® El + E. 



3 
2 



El g) Ef + El g) E^ - 



El 



Ei + El 



E' + Ei 



E', 



'El 



The relations between the z-variable, momentum and the rescaled rapidity u 
transform in the hmit g into 



p = 2z , u = cot z = cot 



P 



(9.2) 
gu 

(9.3) 



Surprisingly enough, this S-matrix cannot be written in the difference form, i.e. as 
a function of one variable being the difference of a properly introduced spectral 
parameter. By construction, this S-matrix satisfies the usual Yang-Baxter equation 



323(^2, Z3)Si3{zi, Z3)Si2{zi, Z2) — Si2{Zi, Z2)Si'i{Zi, Zj,)S23{z2, Z^) , 



(9.4) 



as one can also verify by direct calculation. On the other hand, at one-loop there 
is another "canonical" S-matrix which is a linear combination of the graded identity 
and the usual permutation: 



Qcan 
^12 



Ui - U2 



•-'12 



2i 



-P^ 



12 



(9.5) 



Ui — U2 ~ 2i Ml — ii2 — 2i 

This S-matrix satisfies the same Yang-Baxter equation (9.4). 

The results of [47] imply that the two one-loop S-matrices, (9.2) and (9.5) are 
related through the following transformation 

5"=^°(zi, z,) = U2{z,) [1^(^1)1^2(^2)^12(^1, Z2)Vf\z,)V,-\z2)] U{\Z2) , (9.6) 
where we have introduced the diagonal matrices 



f/(z) =diag(e^^e^^l,l), 

V{z) = diag(e*3, e*3, e"*2, e"'2) . 



(9.7) 
(9.8) 
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The transformation by V is just a gauge transformation which always preserves the 
Yang-Baxter equation. On the other hand, transformation by ?7 is a twist, that 
generically transforms the usual Yang-Baxter equation into the twisted one and vice 
versa [47]. Indeed, S^f^ is nothing else but the one-loop limit of the spin chain S- 
matrix [9]; the latter obeys the twisted Yang-Baxter equation [47]. Note also that 
the twist U does not belong to the symmetry group SU(2) x SU(2) of the string 
S-matrix. 

To understand why at one loop the Yang-Baxter equation is preserved under 
the twisting, we first write the Yang-Baxter equation for S^^'^ by using the relation^^ 
(9.5) 

Us{z2)S23U^\zs)Us{z,)S,,U{\z3)U2{Zi)S,2Ur\z2) = 

= U2{zi)SuUr\z2)U3{zi)SisU{\zs)Us{z2)S23U^\zs) , (9.9) 
which can be reshuffled as follows 

U3{z2)S23U2{z,)Us{zi)SisUr\z3)U2\zs)SuUl{z2) = 

= U2{z,)U3{z,)SuU{\z2)Si3Us{z2)S23U{\z,)U^\z3) . (9.10) 

It is clear now that we will get the usual Yang-Baxter equation for S provided it 
obeys the following relation 

[5,[/®[/] = 0, (9.11) 

where U is an arbitrary diagonal matrix. One can easily verify that both S-matrices, 
(9.2) and (9.5), do indeed satisfy this relation. At higher orders in g the relation 
(9.11) does not hold anymore. The corresponding "all-loop" S-matrix (3.24) satisfies 
only a weaker condition 

[S,G®G] = 0, G G SU(2) X SU(2) , (9.12) 

which is nothing else but the invariance condition for the string S-matrix. As a 
consequence, the Yang-Baxter equation is preserved by the twist transformation 
only at the one-loop order. 

As a final remark, we note that it would be interesting to understand how the 
derivation of the Hirota difference equations for the canonical S-matrix [72] could be 
extended to the "twisted" S-matrix (9.2). This might shed some light on construction 
the fusion procedure for the all-loop S-matrix (3.24). 

^^The gauge transformation by the matrix V decouples from the Yang-Baxter equation. 
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9.3 BAE with nonperiodic fermions 

9.3.1 Bethe wave function and the periodicity conditions 

In any asymptotic domain Q with xq^ <^ xq^ <^ ■ ■ ■ <^ xqj^ where = and 
Qi, . . . , Qtv is a permutation of 1, 2, . . . , A^, the wave function of particles with 
flavors 11,12-, ... can be written as a superposition of plane waves with momenta 
Pi> P2> ■ ■ ■ > Pn 

^t.,^{x,, ...,x^) = J2 ^n'l e^'"-^'^ , (9-13) 

V 

where the sum runs over all permutations of the momenta pi. The scalar product 
Pp ■ xq is deflned as p-p ■ xq = J2k=i Pvf.-'^Qf. ' ^^'^ permutations V and Q 

it satisfies p-p- xq = Ppq-i ■ xx = J2k=i P, o-ii ^'^ where X is the trivial permutation. 

The amplitude ^^'^j^^ is related to the probability of finding the particle with 
the fiavor (the z^-th particle in what follows) carrying the momentum p^^^_-^^ at 
the position Xk. That means that the index ik is attached to the coordinate Xk. As 
a result the wave function (9.13) should satisfy the following symmetry condition for 
any two indices fc, m 

^g...,,...i,„...i^(a;i, Xfc, x™, xn) = 

= (-)^»'=^-^'i;*.-S.....^(a;i.---,^-,---,^A:,...,XAr), (9.14) 

where Vkm is the permutation of k and m, and = if the z-th particle is boson 
and ej = 1 if the i-th particle is fermion, that is one takes the minus sign if both the 
ik-th and Zm-th particles are fermions, and the plus sign otherwise. 

In any two domains Q and Q the amplitudes v4^'.^j^ and -4,^'.®^^ of the same 
plane wave (that is pp -xq = Pp- Xq) are related through the S-matrix. The relation 
can be easily found by representing the amplitudes as the following products of the 
ZF operators 

^X-±<M)---^aM.)^ (9-15) 

and then by using the ZF algebra to relate the amplitudes in the domains Q and 
Q. The +/— sign in this formula is related to the even/odd number of permutations 
of fermions by the permutation Q. To understand the origin of this formula let us 
recall that the indices are attached to the coordinates Xk which explains the order 
of A\ . The dependence of A] of the momentum follows from the coupling Pp^xq^ 
in the exponential of the wave function (9.13). 

To proceed it is convenient to use matrix notations. We introduce the simple 
permutation P12 = Ej which permutes the spaces Vi and V2 but does not 

touch the momenta pi so that S21 = Pi2S{p2,Pi)Pi2 , the graded permutation = 
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{—ly^'^^Ej^E^-, and the graded two-particle S-matrix S'fg which can be written in the 
form = -^i2'S'i2 where = {—ly^^^El (g) Ej is the graded identity. We also define 
S21 = Pi2S{p2,Pi)Pi2li2 = Pi2S{p2,Pi)Pi2 SO that the unitarity condition 8(2821 = I 
is fulfilled. 

Then we multiply the wave function (9.13) and (9.15) by the tensor product of 
N rows ® ® . . . ® E'^ = [E^E"^ ■ ■ ■ e^Y''^'""'' , and (9.13) takes the form 



V 



where 



and the index Qk refers to the location of the row E^'', and Iq is the product of 
graded identities which can be found by representing the permutation Q as a product 
of Y simple permutations Vkm- Q = Vkrmi ■ ■ ■ VkymY, and then J| = Jf^^^ • • • /f^^^. 

Now the ZF algebra can be used to express the amplitudes A^^^ with Qq = V'^Q 
fixed in terms of the amplitude A'^^^°. In particular the amplitudes ^■^'■^ are ex- 
pressed in terms of the incoming amplitude A'^^^ ~ ylj(pi) ■ ■ ■ A^j^{pn). The corre- 
sponding terms in the wave function can be used to derive the periodicity conditions. 

To find the relations, it is convenient to represent 

^ A\{p,) ■ ■ ■ A^M ■ (2o)i...Tv/L = ■ ■ ■ AUpN)iV-'Q),...^pQ^ , 

where {QV^^)i...n is the permutation matrix that acting on the tensor product E''^^ ® 
■ ■ ■ ® E^'^ produces E^'^ ® ■ ■ ■ ® E^'^ . Now we use the ZF algebra to find the relation 

^^is = 4 ■ ■ ■ Aj, ■ 8r,...rAp,^,. . . ,p,J{QV-%..mI'q = 

^xiQo i9^^^^q-^v)i...n8v,...vAPv,,- ■ ■ ,P^^){QV-%..nI'q , 
where S'pj...p^(pp^, . . . ,Pp^) is the multi-particle S-matrix. 
In particular, we find that 

where 8!p_^___p^{p^^, . . . ,p^^) is the graded multi-particle S-matrix. Note that it is 
not a product of two-particle graded S-matrices. 

This formula can be used to find the set of periodicity conditions. We write the 
part of the wave function with the plane wave 6*^'=^* 



V 



A^\^ 5|>,...P,K' • • • e^"---e(xp, < . . . < XV,) . (9.17) 



V 
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The periodicity conditions read 

^(xi, . . . ,Xk = 0, . . . , xn) = ^(xi, . . . ,Xk = L, . . . , XN)Wk , 

where the diagonal matrix W is equal to the identity matrix if the fermions are 
periodic, and it is = {—ly^El if the fermions are anti-periodic. For the su(2|2) 
case we have W = = diag(l, 1,-1,-1) for anti-periodic fermions. 

By using eq.(9.17), we get 
^(xi, . . . = L, . . .,xn) = 

^rp,Lj^x\x J2 S^^...^^_Jp,^,...,p,^_^,pk)Wke^^--^-e{xr, < ... <Xp,_J, 

V:VN=Pk 

Comparing the terms, we obtain 

^^1^ . . . - e^^'=^^|,^...^^,(p,^, . . .,p,^,pk)Wk) = . (9.18) 

To compute the S-matrices, we use their definitions 

(Pv^ ) ■ ■ ■ ^Pvr, ) = A---^N- SkV2-VM (Pk, , • • • , ) , 

(Pv^ ) ■ ■ ■ ^Pv^ )AiPk) = A---^N- Sv2-VMk{Pv2 ,...,p^^,pk) . 

Then we use the ZF algebra to order the product A^,^ {p^^ ) ■ ■ ■ (p^^ ) 
and finally we get the multi-particle S-matrices 

^i^V2^Pv2) ' ' ' ^PAr(Ppjv) = ^1 ■ ■ ■ ■ Sk,k-lSk,k-2 ■■■Ski - Sp^...pj^ 

A^2^Pv2) ' ' ' ^P]v(^'l'jv)^I = ^1 ■ ■ ■ ■ Sk+l,kSk+2,k ■ ■ ■ SNk ■ S-P2-VN ■ 

Thus, for S-p^...-pj^ = 1 eq.(9.18) takes the form 

A^\^ {Sk,k-i ■ ■ ■ SkJlk-i ■ ■ ■ III - e''''^Sk^i,k ■ ■ ■ SNklUk ■ ■ ■ iLWk) = (9.19) 
or, equivalently, 

^x|x (g.p,L _ . . . . . . js^WkHr, ■ ■ ■ Ilk+i^kN ■ ■ ■ Sk,k+i) = . (9.20) 

It is possible to show that the same equations follow if Sp^...pj^ ^ 1 which uses the 
identity Skmlknlmn — IknImnSkm, and also that the terms in the wave function with 
the plane wave e*^^'^^ lead to the same equations. 
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The consistency condition for the system of equations (9.20) requires that the 
matrices 



Tk = Sk,k-1 ■ ■ ■ Sklll^k-l ■ ■ ■ Hl^kl^N ■ ■ ■ Ik,k+l^kN ■ ■ ■ Sk,k+1 

mutually commute. Naturally, we expect that the matrices Tk should be related to 
the monodromy matrix 

T{pa) = -StTAWAS{]^{pA,PN)S{]sf_iipA,PN-l)---S{^{pA,Pl), (9.21) 

where 5*^^ is the fermionic i?-operator defined, e.g., in eq.(102) of [73]. The authors 
of [73] use index notations to define the operator. It is more convenient, however, to 
use the matrix notations and the usual convention for Sjk to work with the operator. 
One can check that it can be written in the following form 

of („. r,,) - i ^j-N^k-N IjkSjk{Pj,Pk) I-.-nII-n if J < ^; /q oo] 

y^j-N^k-N "^jkKPj.Pkj^jk^j-N^k-N 11 J > ■ 

Here /^^ is the graded identity and 



t9 rg j9 ,,79 



To prove the formula, one should use the following representation for the graded 



projection operators E'^^ eq.(28) of [73] 



_ t9 tpP t9 

There are two natural choice for the index A in (9.21), that isy4 = 0oryl = A^ + l. 
The choice leading to appears to be A = + 1 > A;. Then we get 



N ■ 



SikiPA^Pk) = Ik:.NSAk{PA,Pk)lAkH- 

Now we compute the following product 

S{kiPA,Pk)S{^k-liPA,Pk-l) = 

Ik-N^Ak 1% ^I-nH-I-N SA,k~l ^A,k~l ^k-l-N = 

tQ jQ jg q t9 t9 q j9 t9 t9 

^k--N^k-l--N^k~l,k '^Ak J- Ak^ k-l,k'^ A,k~l J- A,k-1 ^k-N ^k-l-N ~ 

^k--N-^k-l--N-^k-l,k S AkS A,k-1 I^k ^A,k-1 ^k~l,k H--N H~h--N ' (9.23) 

where we used the identity 

SA,k-l{pA,Pk) lLl,k ^Ak = H-l,k ^Ak SA,k-l{pA,Pk) ■ 
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The following generalization of the formula (9.23) can be proven by using the math- 
ematical induction 



Sik{pA,Pk)S{ f,_^{pA,Pk-l) ■ ■ ■ S{ f^_JpA,Pk-n) = 

~ ^k---N ^k-n--N ^k-lMk-2---k ^k-n---k^ W-^^J 
X SAk ■ ■ ■ SA,k-n ^Ak ' ' ' ^A,k-n ' ' ' ^k-n---N ^k-l,k^k-2--k ' ' ' H-n--k ■ 

To get the monodromy matrix, we set A; = and n = — 1 in this formula, and 
using the identity 

tq jg jg jg jg j 

J-N-1--N ' ' ' ^N-l,N'^N-2-N ' ' ' ~ ^ ' 

we find the following drastic simplification 

T{pa) = -StrA Wa San ■ ■ ■ Sai Ian ' ' ' • 

Now we choose pa = Pk and use the fact that SAk{Pk,Pk) = —PAk- Recalling that 
our goal is to show that Tip^) = T^, we have 

TiPk) = Str^ Wa San ■ ■ ■ SA,k+i PAk SA,k-i ■ ■ ■ Sai Ian ' ' ' ^ai 

= StiA PAk Wk SkN ■ ■ ■ Sk,k+i ■ SA,k~i ■ ■ ■ Sai Ian ' ' ' Iai 

= Sir A PAk SA,k-i ■ ■ ■ Sai -^A,fc-i ' ' ' ^ai ' Sun ■ ■ ■ Sk,k+i ■ Ian ' ' ' ^Ak ■ 

Now we use that 

iJkN ■ ■ ■ i~>k,k+l ■ ^AN ' ' ' ^Ak ~ ^ AN ' ' ' ^ Ak ' '-'^N ' ' ' i~>k,k+l 

to get 

TiPk) = Str^ Sk,k~i ■■■Ski Il,k~i ■ ■ ■ Hi^kN ■ ■ ■ H,k+i ^AklAk SkN ■ ■ ■ Sk,k+i ■ 
The supertrace can be easily taken 

StrA PAkl'Ak = Tr2 ((-1)^^/ ® i?:) (K ® K) {{-ly^'^Ej ® E^) 
= (-1^+^^^: = /, 

and, therefore, we show that T{pk) = Tk- Since T(u)T{v) = T{v)T{u) for any u and 
V, we have shown that the periodicity equations (9.20) are consistent. 

^^In framework of the algebraic Bethe Ansatz twisted boundary conditions for Hubbard-like 
models have been studied in [631. 
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9.3.2 Two-particle Bethe equations 

To see how the formulas of the previous subsection work let us consider a two-particle 
wave function given by 

1^2121 • , • 2^1121 •, • . (9.25) 

According to (9.15), we can identify 

It is clear that the amplitudes Alf^"^ and A^j^"^^ correspond to the in- and out-states, 
respectively. By using the ZF algebra we find 

A]iP2)4iPi) = S'Ji{p2,Pi)AI{pMKp2) 4'"' = (-)"^^^?(P2,Pi)4?"' • 

In a similar way we get 

4f ' - 4iP2)A]ipi) , - i-r^A]{p,)4ip,) , 

and 

A]{pi)Ai{p,) = 5j^(pi,P2)4(P2)4(Pi) = i-r'^si^{p,,pMT ■ 

12 1 12 21 1 12 

The amplitudes A^^j and A^j are not independent. By the symmetry condition 
(9.14) they are related to each other as follows 

^J'" = (-)-^?" = s-{p.,P2)^'^ =, A^ = s^:;{p,^pmT ■ 

The wave function (9.25) can be written in the matrix form by multiplying it by the 
row -E* ® and summing over z, j. Then we get 

(A (e^Pi"i+*f2X2 + S^^p^^ ^ip2X^+^p^X2^ < 
[XuX2) <^J^^p9^^ip,x2+^P2X^ ^s9^^ip2X2+^p^x^^ X2 < Xi ' ^ ' 

where we recall that Pfg = {—ly^'^^Ej ® is the graded permutation and S21 is the 
graded S-matrix S^^ = Pi2S{p2,Pi)Pi2lf2 = Pi2S{p2,Pi)Pf2- 
The (quasi)-periodicity condition can be easily imposed 

^(Xi,X2) = ^(Xi +L,X2)Vri, '^{Xi,X2) = 'if{Xi,X2- L)W2, Xi < X2 , 

where the matrix W is equal to / for periodic boundary conditions and to S for 
anti-periodic boundary conditions for fermions. By using the wave e*^*^'' this leads 
to the following equations 

^ (1 - e'P'^Si^Wi) = , ^ (1 - e-'P'^Sl^W2) = , 
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or by using the wave e*^2a;2+«pia;2 

A {S21P12 - e'P^'^Pf^Wi) = , A {S21P12 - e-'P^'^Pf^W^) = . 

These two sets of the periodicity conditions are obviously equivalent because Pf2^i = 
W2Pf2- Let us also mention that the equations are compatible if the matrices WiSf2 
and S21W2 commute, and this follows from unitarity SI2S21 = ^ and 

W1W2SI2 = SI2W1W2 . 

Let us now see how the nesting procedure works for the case of one A\ boson 
and one A\ fermion. Consider the system of equations 



A 



13 ~ 

AS^"" = SlUP2,Pi)AlT' + Sllip2,Pi)A,, 



13 ' 
12|12 



(9.27) 



Assuming that S^j are matrix elements of the string S-matrix S, we get 



^2i|i2 ^ So{p2,Pi] 
AI?^ = So{p2,pi] 



^1 -^2 c?Pi^/^12| ^ 

xt - X2 v{Pi) eiP^ 



2\i2 xt -x^ r]{p2)eiP' 12112 



Xt-X2 ^(Pl)^12|12 , X'^-x'[^_ 



(9.28) 



x+ ri{p2 



31 



+ 



12|12 
13 



where 5*0(^1,^2) is the scalar prefactor. 

For the amplitudes of interest the general Bethe equations 

g-ipiL^12|12 ^ 



-1 



SfAP2.P1) A 



12|12 
kl 



(9.29) 



read as follows 



12|12 _ ^21|1 



^_.,,,^i2|i2 ^ sII{p2.pMT + S^}{p2.pMn 

e--M-i- = {-ir\sll{p2,pMT + S!1{P2,pMT 



2 

31 ' 



21|12 
13 ' 



(9.30) 



where we have used eqs.(9.27). Note that in eq.(9.29) the multiplier (—1)^^* takes 
into account the boundary conditions for fermions: e = for periodic fermions and 
e = 1 for anti-periodic ones, respectively. 

The system (9.28) can be solved in two different ways depending on the choice 
of the first level vacuum [9]. Below we present both solutions. 

• Regarding Ai . . . Ai as the first level vacuum, we first choose the following 

ansatz 



aT = fiP2)S{pi) , 



A'T = Sllip2,p,)f{p^)S{p2) 



Air = fiPi) , 

Afr = Sllip2,P,)fip2), 



(9.31) 
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where Sll{pi,p2) is the corresponding element of the string S-matrix. One can 
easily show that this ansatz indeed solves the system (9.28) provided we take 

, , e'2 x+ -X- i^y-x- 
/(P) = S{p) = e2 



vip) y ~ X ' 2/ — x+ 

According to eqs.(9.31), the last formulae give 

e-*^^V(P2)5(pi)=5i\i(l>2,Pi)/(P2), 

e-^^^^/(pi) = i-lrSllip2.p^)fip^)S{p2) 

and we derive the corresponding Bethe equations 

e'P^'^ = Sll{p„p2)S{p,), 
i-ir = Sip,)S{p,). 

If we choose A3 ... A3 as the first level vacuum, we modify the ansatz for the 
corresponding amplitudes as follows 



AT = fiPi) , AT = fiP2)S{pr) . 

AT = S!Up2,Pi)fiP2) , AT' = S||ip,,p^)fip^)Sip,) . 

Note that S'||(pi,p2 ) = —So{pi,p2)- This time satisfaction of eqs. (9.28) requires 
one to choose 

fip) = ^(P)e-^'^^ , S{p) = -e-i^^ . 

y — x~ y — x~ 

The Bethe equations (9.30) read 

e-'^^^/(pi) = ^l3(P2,Pi)/(pi)^(p2), 
e-^^^^/(p2)5(pi) = (-l)^5i(p2,Pi)/(P2), 

and, therefore, we find 

e^P^^ = {-iysll{p^,p,)S{p,) = (-iySo{pi,P2)^l^e~^'^ 



2 

xi -y^ 



[-ly = s{p^)s{p2) . 



This completes consideration of our simple example illustrating the dependence of 
the Bethe equations on the periodicity conditions for fermions. 



9.4 Large/small g expansions of solutions to the bound state equation 

The four general solutions of the bound state equation (7.8) are 

(y/f?^ sin^l + l + lj [cos | ,J g'' sin^ f + 1 ± y^cos^ f - g"^ sin^ l) 

= ' ^^-^^^ 

iJg'' sin^l + i - i) fcos f ./Tsh^^iTI ± ^/^TT^^^^T) 

e-^ = ^ \.2v ' (9-34) 

g^ sm I 

where only the first two solutions (9.33) correspond to states with positive energy. 

The large g dependence of q of the bound state solutions with momentum ex- 
ceeding Per is obtained by expanding (9.33) in powers of with the bound state 
momentum p kept fixed^^ 

" ~ %W| ^ ' U " 2^?2 sin3 I ) + ^ • ^^-^^^ 

To find the large g dependence of q of the bound state solutions with momentum 
smaller than one should take into account that '^1 \f9 as oo, and 

therefore one should consider a bound state with momentum p of the order 1 / 
and keep the product fixed in the large g expansion 



l±A/l-2jf- 4 / 



?± = 2 - 1 - ^ ± • (9.36) 



gp 'ig^y^ \ 16 



1 _ pV 

16 



The small g dependence of q of the bound state solutions with momentum smaller 
than Per is obtained by expanding (9.33) at small g with the bound state momentum 
p kept fixed 

4 cos ~ T) 

g+ = -2 log (7 + log — ^ + ^^(1 + 3cosp) tan^ ^ + 0(/) , (9.37) 

Sill 2 

2 

g_ = - log cos I + ^ sin^ | tan^ | + 0((74) . (9.38) 

To find the small g dependence of q of the bound state solutions with the mo- 
mentum exceeding pcr, one should take into account that pcr tt — 2g as g ^ 0. 
Then, one can parametrize p as follows 

p = TT — 2g cos a , 

and keep a fixed in the expansion. Then we get 

= -log^ + — (2 + cos2a) ±i (a + — sin 2a - — cot + O(g^) . (9.39) 
2 4 V 6 6 / 



25 



We assume here that p e (0,7r). 
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